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1.  Electric  Vacuum  Gyroscope 

The  analysis  of  slow-down  and  precession  torques  resulting  from  radial  mass 
unbalance  of  an  electric  vacuum  gyroscope  rotor  Is  described.  Continued  exper¬ 
iments  on  the  verification  of  torques  resulting  from  deformed  rotor  translational 
misalignment  are  described  and  the  results  given. 

2.  Surface  and  Atomic  Physics 

The  pressure  shift  of  the  hyperfine  structure  of  atomic  nitrogen  N14  in  the 
ground  state  has  been  Investigated  for  the  following  buffer  gases:  molecular 
nitrogen,  helium,  neon  and  argon. 

3 .  Computers 

Operation  statistics,  logic  modifications,  and  programming  efforts  are  re¬ 
ported  for  the  CSX-1.  For  the  CDC  l6o4,  operation  statistics  are  given.  Studies 
in  which  a  computer  plays  both  roles  in  a  game  are  reported.  Provisions  by  which 
a  computer  may  produce  multivoiced  music  is  reported.  An  abstract  of  a  recent 
report  on  Automatic  Air  Traffic  Control  is  presented. 

4 .  System  Uteory 

Various  results  related  to  control  optimization  and  feedback  theory  were 
obtained.  Various  computer  programs  necessary  for  contact  network  synthesis, 
self-diagnosis,  and  graph  theory  were  completed.  Further  results  on  self -repair 
and  on  communication  nets  were  also  obtained. 

3.  PLATO 

In  January  1963,  ILLIAC  was  removed  from  service.  The  computer  control  for 
PLATO  was  switched  from  the  ILLIAC  to  the  CDC  l6o4  located  at  the  Coordinated 
Science  Laboratory.  Two  teaching  logics  using  PLATO  Laboratory  were  completed. 

One  is  used  for  demonstration  purposes  to  illustrate  the  teaching  of  Archimedes 
Principle  by  simulated  laboratory  experiments.  The  other  logic  will  be  used  for 
inquiry  training.  Preliminary  analysis  of  student  electrocardiogram  data  indicates 
that  student -discovery  is  reflected  in  the  heart  rate.  Testing  of  the  PLATO  III 
equipment  was  started  this  quarter. 


6.  Vacuum  Instrumentation 
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Study  of  sticking  probability  of  Ng  using  partial  pressure  gauges  has  con¬ 
tinued.  A  study  of  the  anomalous  behavior  of  the  Bayard -Alpert  gauge  at  very  low 
pressure  vas  begun  this  quarter,  using  both  partial  and  total  pressure  gauges. 

7 .  Plasma  Physics 

Although  Nordsieck's  method  of  solution  of  the  Boltzmann  equation  for  a 
shock  wave  appears,  by  most  criteria,  to  be  successful,  the  Monte  Carlo  calcula¬ 
tion  of  the  collision  integral  is  biased.  A  series  of  new  routines  has  therefore 
been  developed  that  will  permit  efficient  study  of  the  collision  integral  calcu¬ 
lations.  The  momentum  transfer  cross  section  for  electrons  colliding  with  Neon 
atoms  was  measured  as  a  function  of  temperature  in  the  range  200  -  600°  K.  The 
results  indicate  a  minimum  in  the  cross  section  at  0.11  eV.  An  experiment  has 
been  performed  to  demonstrate  the  incoherent  scattering  of  microwaves  from  a 
plasma.  A  resonance  was  found  when  the  frequency  of  the  scattered  radiation 
corresponded  to  the  sum  of  incident  frequency  plus  plasma  frequency.  A  new  im¬ 
proved  model  of  the  plasma  delay  line  has  been  constructed  that  should  result  in 
a  better  focusing  of  the  electrons  and  provide  more  experimental  flexibility* 

8.  High  Magnetic  Field  Superconductors 

Measurements  have  been  made  of  the  critical  current  as  a  function  of  temper¬ 
ature  and  magnetic  field  of  single  crystalline  vanadium  and  compared  with  measure 
ments  on  similar  material  which  has  been  cold  worked.  Clusters  of  Nb^Sn  crystals 
have  been  grown  from  saturated  solutions  of  Nb  in  Sn  at  temperatures  between 
1,000°C  and  1,250°C.  Some  developments  in  fabrication  of  Nb^Sn  wire  have  been 
made  in  cooperation  with  the  Fansteel  Metallurgical  Corporation. 

9.  Vacuum  Electrical  Breakdown 

The  preliminary  measurements  previously  reported  have  been  repeated  with 
greater  care.  The  electron  microscope  studies  have  been  continued. 

10 .  Physical  Electronics 

The  voltage -current  characteristics  of  metal -oxide -metal  sandwiches  have 


been  studied 
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11,  Arms  Control  and  Disarmament  Beaearch 

Beginning  part-time  studies,  In  cooperation  with  faculty  from  other  depart 
ments,  are  reported. 
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1.  ELECTRIC  VACUUM  GYROSCOPE 


W.  Khoebel 

R. 

W.  Anderson 

V. 

Krone 

0.  Skaperdas 

J. 

D.  Gooch 

D. 

Lee 

Allen 

B. 

D.  Kirkwood 

W. 

C.  Prothe 

1.1  Introduction 

The  electric  vacuum  gyroscope,  a  high-precision  inertial  navigation  instru¬ 
ment  employing  a  coasting  spherical  metal  rotor  supported  without  physical  contact 
in  vacuum  by  electric  fields,  is  being  investigated.  The  objectives  are  (l)  to 
obtain  a  complete  understanding  of  subtle  effects  which  may  influence  the  design 
of  ultra -high  precision  electric  vacuum  gyroscopes  and  (2)  to  develop  new  tech¬ 
niques  for  gyro  improvement. 

Long  term  experiments  using  two  different  gyro  test  setups  are  being  conducted 
to  experimentally  verify  predicted  small  changes  in  performance  resulting  from 
various  parameter  adjustments. 

1.2  Gyro  Torque  Analysis 

A  complete  analysis  of  the  sources  of  precession  and  slow-down  torques  on  the 
electrically  supported  gyro  should  include  the  effects  of  rotor  mass  unbalance1. 

The  mass  center  and  geometric  center  do  not  coincide  in  general,  and  a  vector 
specifying  the  location  of  the  geometric  center  with  respect  to  the  mass  center 
has  two  components,  one  parallel  to  the  spin  axis  and  the  other  normal  to  it.  Con¬ 
sideration  of  the  two  is  best  done  separately. 

The  parallel  component  produces  the  classical  mass -unbalance  precession  torque 
which  appears  because  acceleration  forces  on  the  rotor  converge  at  the  mass  center 
while  the  electric  forces  converge  at  the  geometric  center,  and  a  torque  couple 
is  created.  The  analysis  of  the  effect  of  this  torque  on  gyro  performance  is 
straightforward  and  is  not  discussed  here,  but  the  effect  of  the  radial  component 
is  more  subtle. 

Precession  torques  may  occur  indirectly  because  of  radial  unbalance.  The 
rotor  tends  to  spin  about  an  axis  which  passes  through  the  mass  center,  but  the 
support  servos  tend  to  force  the  rotor  to  spin  about  an  axis  through  Its  geometric 

1  Minneapolis -Honeywell  Reports:  1515-QR1  (July,  1959),  1901-QR4  (March,  i960), 
1901-TR1  (April,  1961). 
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center.  In  general,  the  spin  axis  passes  through  some  other  point  which 

varies  with  the  spin  speed:  at  very  low  speeds  the  servos  win  control,  and 
the  rotor  spin  axis  passes  through  the  geometric,  and,  at  very  high  speeds, 
the  servos  lose  control  and  the  spin  axis  passes  through  the  mass  center. 

First,  assume  an  infinite  number  of  electrodes  and  let  each  diametrically 
opposing  pair  have  Its  own  push-pull  support  servo.  As  the  spherical  rotor 
spins  about  some  axis  not  through  the  geometric  center,  each  electrode  pair 
"sees"  a  sinusoidally  varying  gap.  Most  servos  sense  the  gap  through  the  accom¬ 
panying  capacitance  of  the  electrode  to  the  rotor.  The  sensed  signal  Is 
amplified,  frequency  compensated  to  provide  stability,  and  then  converted  to 
a  correction  voltage  which  is  applied  to  the  electrode  which,  in  turn,  applies 
a  corrective  force  on  the  rotor.  The  servo  feedback  frequency  response  Is 
specified  by  a  gain  or  stiffhess  (force  per  unit  gap  change)  and  a  phase  angle 
between  output  force  and  gap  change. 

If  the  frequency  response  of  all  servos  Is  Identical,  then  an  analysis 
shows  that  there  are  no  precession  torques,  but  there  Is  the  slow-down  torque. 
The  Integrated  value  of  this  torque  is 

T  =  -  3  Kp2a  sinQ, 

where  T  is  the  torque  (dyne  cm), 

K  is  the  servo  gain  (force  per  unit  area/fractional  gap  change), 

p  is  the  radial  distance  between  actual  spin  axis  and  a  line  through 
the  geometric  center  (cm), 

a  is  the  sphere  radius  (cm),  and 

C  Is  the  servo  phase  angle. 

At  rotor  speeds  where  the  servo  output  response  leads  the  gap  change  (Q  >  0) 
the  rotor  slows,  but  when  C  <  0,  the  rotor  accelerates.  It  follows  that  If  all 
support  servos  are  not  Identical,  precession  torques  appear. 

Consider  the  more  practical  case  of  electrodes  placed  symmetrically  over 
the  sphere,  either  In  the  pattern  of  the  regular  polyhedron  or.  If  a  large 
number  of  electrodes  is  desired,  In  the  pattern  obtained  when  the  spherical 
faces  of  the  regular  polyhedrons  are  broken  up  Into  sub-faces.  Consider  the 
torque  from  one  electrode,  which  may  be  written  as  two  torques— a  slow-down 
and  a  precession  torque: 


3- 


T.lo»  tom  -  -  5  k"2“ 

^.CUlOD  ■  '  E  Kp2* 

where  Q  is  a  function  of  the  solid  angle  subtended  by  the  electrode, 

Y  is  the  angle  of  the  spin  axis  from  a  radial  line  to  the  center  of 
the  electrode, 

i  is  a  unit  vector  in  the  plane  in  which  Y  is  measured  and  normal 
~  to  the  spin  axis, 

and  is  a  unit  vector  normal  to  the  plane  in  which  Y  is  measured. 

Both  torques  are  functions  of  the  angular  location  of  the  spin  axis  with 
respect  to  the  electrode,  and  both  functions  are  symmetric;  i.e.,  the  torques 
from  one  electrode  and  the  electrode  diametrically  opposite  add.  After  adding 
the  contribution  of  the  slow-down  torque  from  all  electrode  pairs,  the  total 
torque  is  less  than,  but  approaches,  the  torque  calculated  in  the  infinitely 
many  electrode  case. 

Die  individual  electrode  precession  torques  are  of  the  same  form  as  the 
second-order  rotor  distortion  torques,  except  that,  in  this  case,  their  magni¬ 
tudes  and  directions  depend  on  the  servo  gain  and  phase  parameters,  which  are 
a  function  of  spin  speed.  The  total  torque  will  vanish  if  the  frequency 
characteristics  of  all  support  servos  are  identical. 

A  complete  analysis  of  this  work  is  being  included  in  a  comprehensive 
report  on  electric  vacuum  gyroscope  precession  torques  by  D.  Allen  to  be  pub¬ 
lished  during  the  next  quarter. 

D.  Allen 

1.3  Gyro  Rotors 

It  was  reported  in  the  last  quarter  that  degassing  from  two  rotors  man¬ 
ufactured  by  the  American  Beryllium  Company  prevented  their  use  in  a  gyro 
system •  Mild  bakeouts  did  not  clean  up  the  rotors,  nor  did  sonic  cleaning 
of  one  of  them  decrease  its  degassing  rate. 

Two  other  beryllium  rotors,  also  fabricated  by  American  Beryllium,  were 
sent  to  Speedring  Corporation  of  Warren,  Michigan,  for  balancing.  After  a 
brief  attempt  at  balancing,  these  two  rotors  were  returned  to  us  with  a 
notation  that  many  small  surface  pit  holes  were  being  uncovered  during  lapping 
and  that  these  pits  loaded  with  lapping  compound  and  caused  a  tearing  condition. 
The  previously  mentioned  rotor  degassing  may  be  due  to  minute  pit  holes  which 
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have  trapped  lapping  compounds  or  other  machine  shop  residues.  Regardless 
of  its  nature,  this  is  not  a  basic  problem  since  excellent  rotors  are 
obtained  from  Minneapolis -Honeywell. 

1.4  Gyro  Test  Data 

Previous  runs  concentrated  on  measuring  the  effects  of  a  constant  dis¬ 
placement  of  a  centrifugally  deformed  rotor  -with  respect  to  its  electrodes 

2 

in  the  ji  and  -jL  direction  •  It  was  shown  that  the  experimental  results 
corroborated  the  theoretical  predictions;  namely,  that  for  this  type  of 
rotor  displacement  the  spin  axis  locus  should  show  spiraling  and  culmination 
shifts,  with  the  closure  time  not  affected.  The  theoretical  results  were 
shown  in  graphical  form  in  Fig.  1  of  the  previous  progress  report.  As  a 
summary  for  the  series  of  runs  56b,  C,  and  D,  the  predicted  ratio  of  the  cul¬ 
mination  shift  to  the  spiraling  rate  per  day  for  a  15°  coning  radius  about 
north,  as  obtained  from  Fig.  1  of  the  previous  progress  report,  is  420/350  = 

1.2*  The  ratio  actually  measured  was  1.3  +  0.1,  the  error  being  attributed 
to  the  difficulty  in  accurately  determining  the  culmination  point.  The  spir¬ 
aling  and  culmination  shift  directions  also  agreed  with  theory. 

In  the  past  quarter,  an  excellent  beryllium  rotor,  fabricated  and  bal¬ 
anced  by  Minneapolis -Honeywell,  was  tested.  In  the  series  of  runs  63  A,  B, 
and  C  the  centrifugally  distorted  rotor  was  given  a  constant  displacement 
with  respect  to  its  electrodes  in  the  ^  and  directions  (see  Fig.  1,  Prog¬ 
ress  Report  for  Dec.,  1961,  Jan.  and  Feb.,  1962).  The  theoretical  results 
for  this  type  of  displacement,  obtained  by  solving  equations  1.1  and  1.2  of 
the  previous  progress  report  on  a  digital  computer,  are  plotted  in  Fig.  1.1. 

It  is  to  be  noted  here  that  the  closure  time  and  elongation  shifts  are  affected 
and  that  no  spiraling  takes  place.  The  curves  show  the  predicted  elongation 
shift  and  closure  time  difference,  from  one  sidereal  day,  versus  the  nominal 
angle  at  which  the  spin  axis  cones  about  north.  The  curve  for  elongation  shift 
is  antisymmetric,  and  that  for  the  closure  time  difference  is  symmetric  about 
the  90°  point.  In  runs  63  A,  B,  and  C,  the  nominal  coning  radius  was  22.5°. 

From  Fig.  1.1,  the  predicted  ratio  of  the  closure  time  difference  to  the  elong¬ 
ation  shift  for  a  coning  radius  of  22.5  degrees  is  5320/305  =  17*4.  The  actual 
ratio  measured  was  17.2  +  1.2,  the  error  again  being  attributed  to  the  difficulty 


Fig.  1,  Progress  Report  for  Dec.  1961,  Jan.  and  Feb.,  1962 


in  rad/day 
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In  accurately  determining  the  elongation  point.  The  directions  of  the  elonga¬ 
tion  shifts  and  the  closure  time  differences  agreed  with  those  indicated  in 
Fig.  1.1. 

Since  we  have  previously  qualitatively  corroborated  the  effect  of  a  con¬ 
stant  rotor  displacement  in  the  k  direction,  showing  that  it  is  equivalent 
to  a  change  in  mass  unbalance,  we  should  now  be  able  to  predict  the  effects  of 
any  constant  displacement  of  a  centrifugally  deformed  (first  and  second  order 
distortion)  rotor  with  respect  to  its  electrodes. 

The  knowledge  gained  during  the  previous  quarter  was  used  in  an  attempt 
to  "tune  up"  the  beryllium  gyro  rotor  described  above.  In  other  words,  a 
series  of  experiments  were  performed  using  a  theoretically  determined  adjust¬ 
ment  to  correct  for  mass  unbalance,  sag  and  non-centering  of  the  rotor.  The 
rotor  was  successively  given  constant  displacements  with  respect  to  its  elec¬ 
trodes  in  the  i^  and  k  directions  with  expected  changes  in  spiraling  rates 
and  closure  time  differences.  The  rotor  rotation  frequency  was  kept  at  a  low 
110  rps  in  order  to  minimize  damage  in  case  of  a  levitation  failure.  Conse¬ 
quently,  the  drift  rate  due  to  mass  unbalance  was  too  large  to  cancel  out  by 
displacing  the  rotor  a  nominal  amount  along  its  spin  axis  (k  direction).  Changes 
in  closure  time  of  2-l/2  minutes  were  easily  obtained  and  the  spiraling  could 
be  changed  in  either  direction  or  eliminated.  A  failure  in  the  levitation 
system  terminated  this  experiment.  The  effects  of  this  accident  on  the  rotor 
are  now  being  Investigated. 


D.  Skaperdas 
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2.  SUBFACE  AND  ATOMIC  PHYSICS 
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2.1  Pressure  Shift  of  the  Atomic  Nitrogen  Hyperfine  Structure 

The  pressure  shift  of  the  hyperfine  structure  of  atomic  nitrogen  l4  in 
3  4 

the  (2p)  Sj^2  electronic  ground  state  vith  molecular  nitrogen,  helium,  neon, 
and  argon  as  buffer  gases  has  been  Investigated.  Buffer  gas  pressures  ranged 
from  7  Torr  to  &J  Torr.  Die  results  complement  previous  studies  made  both  at 
Harvard"^  and  at  Illinois. 

Oriented  atomic  nitrogen  vas  obtained  through  spin  exchange  vith  optically 
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pumped  cesium  vapor  using  the  flov  technique  described  by  Holloway,  et  al., 
vith  several  improvements.  The  a.c.  magnetic  field  compensation  has  been  im- 

5 

proved  according  to  a  suggestion  from  Lambert,  resulting  in  a  narrower  reson¬ 
ance  line  width,  hence  providing  more  accurate  results.  The  dissociation  of 
molecular  nitrogen  was  effected  by  a  new  gaseous  discharge  utilizing  stainless 
steel  electrodes  powered  by  a  100  watt  1.3  mm/sec  oscillator.  Finally,  a  cesium 
light  source^  with  considerably  greater  D^  intensity  than  previously  available 
was  constructed,  providing  a  much  better  signal -to-nolse  ratio. 

From  all  observed  nitrogen  hyperfine  transitions,  the  least  field  depend¬ 
ent  (F  =  3/2,  mp  =  +3/2  F  =  l/2,  m^,  =  +l/2)  were  chosen  in  order  to  provide 
the  most  accurate  determination  of  the  frequency  shift.  Measurement  were  made 
at  twelve  different  values  of  pressure  for  each  buffer  gas,  with  an  average  of 
ten  readings  per  data  point.  The  average  standard  deviation  for  any  one  of 


R.  H.  Lambert  and  F.  M.  Pipkin,  Bull.  Amer.  Phys.  Soc.  7_,  26  (1962). 

R.  H.  Lambert  and  F.  M.  Pipkin,  Phys.  Rev.  129,  1233  (19^3) • 

W.  W.  Holloway,  Jr.,  E.  Luscher,  and  R.  Novick,  Bull.  Amer.  Phys.  Soc. 

26  (1962). 

W.  W.  Holloway,  Jr.,  E.  Luscher,  and  R.  Novick,  Phys.  Rev.  126,  2109  (1962)* 
R.  H.  Lambert,  private  communication. 

R.  G.  Brewer,  Rev.  Sci.  Instrum.  32,  1356  (1961). 
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these  points  is  less  than  15  cps.  The  results  of  these  neasureaents  are  rep¬ 
resented  In  Table  I,  while  a  comparison  of  the  experimental  results  with  the 

7 

theoretical  values  calculated  by  Adrian  Is  presented  in  Table  II.  All  readings 
with  helium,  neon,  and  argon  as  buffer  gas  were  made  with  the  same  constant 
amount  of  5.0  +  0.5  Torr  of  nitrogen  in  the  system.  The  zero  pressure  inter¬ 
cepts  for  helium  and  neon  must  therefore  be  shifted  down  by  7  cps,  and  6  cps, 
respectively,  while  the  corresponding  slopes  remain  unaffected.  The  corrected 
curves  are  graphed  in  Fig.  2.1.  The  average  of  the  values  for  the  zero  pressure 

splitting  is  15  676  57^  +  10  cps,  compared  with  the  value  given  by  Lambert  and 
2  " 

Pipkin  of  15  676  592  +  6  cps  for  helium  only,  and  the  value  given  by  Holloway, 
et  al.^  of  15  676  590  +  cps,  for  nitrogen  only.  There  is  fairly  good  agreement 
in  this  experimental  data.  The  zero  pressure  splittings  had  not  been  measured 
before  in  the  presence  of  neon  or  argon.  Unfortunately,  we  were  not  able  to 
observe  the  nitrogen  resonances  using  krypton  or  xenon  as  buffer  gas. 

E.  Luscher 

F.  A.  Franz 
D.  W.  Smith 
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F.  J.  Adrian,  Fhys.  Rev.  12£,  837  (1962)- 
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Table  I 


Buffer  Gas 

Pressure  Range 
(mm  Hg) 

vQ(F  =  3/2-F  =  1/2) 
(cycles/sec) 

^dv/dp^> 

(c/sec  mm  Hg) 

He 

19.9  -  86.1 

15,676,380  +  8 

0.18  +  0.056 

Ne 

14. 2  -  87.0 

15,676,379  +  10 

0.365  +  0.060 

Ar 

10.9  -  53-6 

15,676,381  +  8 

2.56  +  0.09 

6.7  -  85.0 

15,676,369  +  10 

2.49  +  0.05 

Values  of  the  zero  pressure  Intercepts  and  slopes  of  the  frequency  vs  buffer  gas 
pressure  curves  for  the  (F  =  3/2  “*  F  =  l/2)  hyperfine  transition,  as  derived 
from  a  least  squares  fit  to  the  experimental  data.  Several  of  the  zero  pressure 
intercepts  must  be  adjusted,  as  is  explained  in  the  text.  Errors  quoted  are 
equal  to  two  standard  deviations  from  the  mean. 

Table  II 

Effect  of  Various  Buffer  Gases  on  the  Nitrogen  Hyperfine  Splitting  Constant 

/  a V*\  /  (dA)  /(M) 

Buffer  Gas  \  dp  /  \  dp  A.  'dp'  '  'dp' 

\  *  /exp  \  F  /theory  *  exp  *  theory 

(c/sec  mm  Hg)  (c/sec  mm  Big) 


He 

0.121 

+  0.037 

1.1 

0.11 

Ne 

0 

to 

$■ 

+  o.o4o 

1.8 

o.i4 

Ar 

1.71 

+  0.06 

5.o 

0.34 

No 

1.66 

+  0.03 

4.8 

0.34 
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5.  COMPUTER 

R.  M.  Brown 

L.  Kypta 
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M.  Levitt 

R.  Trogdon 

J.  Divllblss 

V.  V.  Llchtenberger 

F.  Vise 

J.  Govaia 

P.  Ponzo 

N.  Wax 

R.  Jenks 

S.  Srivastava 

5.1  Introduction 

This  division  is  concerned  with  the  design  and  application  of  digital 
computers  for  information  processing,  primarily  in  the  areas  of  real  time  and 
non -numerical  operations.  The  present  work  lies  in  three  areas:  the  develop¬ 
ment  and  supervision  of  the  CSL  computing  facilities  centering  on  the  CSX-1 
and  a  CDC  l6o4  computer,  cuad  studies  in  artificial  intelligence. 

5 .2  CSX-1  Computer 
5.2.1  Operations 

Period:  November  29,  19^2  to  February  28,  1965 
Total  Running  Time  Vt-9.5  hrs 

Ave.  per  day  (5  day  week)  6.9  hrs 


Scheduled  Maintenance  Time 


0.0  hrs 


Emergency  Maintenance  Time 
Operational  Time 

Breakdown  of  Emergency  Maintenance  Time 


1.2  hrs  (.27  0/0) 
¥»8.5  hrs  (99*73  0/0) 


Cold  Solder  Joint  1  hr 

Replace  6  transistors  .2  hr 

The  computer  was  shut  down  for  five  weeks  during  this  quarter  for  the 
installation  of  the  V  =  7  modification  (see  below)  and  a  temperature  sensing 
unit.  lhe  transistor  failures  mentioned  above  occurred  during  the  installation 
and  checkout  of  the  temperature  sensing  unit.  They  were  in  three  driver  stages. 

5.2.2  Logic  Modifications 

The  V  =  7  modification  affects  those  orders  for  which  the  V  or  n 
bits  refer  to  an  accumulator.  For  these  orders  if  V  or  n  is  Identically  seven, 
the  word  or  words  Immediately  following  the  order  in  the  program  specify 
memory  locations.  The  effect  of  this  modification  is  to  permit  any  memory 
located  to  be  designated  as  an  accumulator  and  processed  accordingly. 

In  addition  to  V  -  7  modifications  the  following  changes  in  the  order 
code  are  being  made. 
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Addition  of  Input  and  output  orders  directly  to  and  from  memory 

Addition  of  tit -compliment  orders 

Removal  of  tbe  two  logical  sum  orders  and  replacing  them  with  two 
special  masked  compare  orders 

Modification  of  the  existing  compare  orders 

Addition  of  a  memory  Increment  order. 

Details  of  these  changes  along  with  the  V  =  7  modification  will  appear  In 
a  new  CSX-1  manual. 

J.  Stifle 
R.  Trogdon 

3.2.3  Programming 

A  new  CSL-1  simulator,  SAC SEC,  was  written  for  the  CDC  l6ok  In 
order  to  accomodate  the  modifications  In  the  order  code.  This  program  enables 
programs  written  In  the  new  order  code  to  be  assembled  and  checked-out  prior 
to  the  completion  of  hardware  modifications. 

Tbe  CSX-1  assembly  program,  SAC,  Is  currently  being  revised  to  make  Input 
format  more  flexible  and  compatible  with  l6ok  assembly  programs .  The  major 
modifications  to  the  assembler,  scheduled  for  completion  during  the  next  quarter, 
permit  the  following: 

1.  The  use  of  symbolic  addresses  consisting  of  up  to  eight 
alphanumeric  characters. 

2.  Mnemonic  codes  of  up  to  nine  alphabetical  characters 
(normally  three ) . 

3.  New  Input-output  "macro-instructions"  similar  to  those  used 
In  Fortran. 

k.  Format  statements  associated  with  these  macro -Instruct Ions 
which  also  may  be  used  for  packing  binary  data  Into  program 
words. 

Routines  to  read  and  write  magnetic  tape,  and  output  to  the  1612  printer 

were  completed.  The  routine  "COPY"  which  takes  Input  from  any  Input  channel 

and  reproduces  It  on  any  output  channel  has  been  completed  and  checked  out. 

R.  M.  Brown 
R.  D.  Jenks 
F.  G.  Vise 
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3. 5  CPC  l6o4  Computer 
?.?.!  Operations 

Period:  December  1,  to  February  28,  1963 

Total  Running  Time:  986.5  hours 

Average  per  day  (five  day  week)  15-91  hours 

Average  per  day  (including  Saturday,  Sunday  and  holidays) 

IO.96  hours 

Operational  Time:  915  hours  (92.55  0/0) 

Preventive  Maintenance  Time:  69  hours  (6.99  0/0) 

Emergency  Maintenance  Time :  4.5  hours  (0.45  0/0) 

Actual  Operating  Time/Scheduled  Operating  Time  99*51  0/0 

Breakdown  of  Emergency  Maintenance  Time : 

Main  frame  of  computer:  1  hour  (solder  short  between  pins 

in  inter-chassis  socket) 

1607  magnetic  tape  unit:  4.5  hours  (tape  wrapped  up  on 

pinch  roller) 

E.  Neff 


3.4  Artificial  Intelligence 

3.4.1  Simple  Game  Playing 

The  preceding  progress  report  described  a  computer  program  which 
enables  the  computer  to  play  a  simple  game  against  a  human  opponent.  This 
program  demonstrated  that  the  machine  could  play  "correctly"  (i.e . ,  according 
to  the  rules)  and,  more  importantly,  could  learn  to  play  a  respectable  game 
against  opponents  of  varying  skill. 

During  this  quarter  a  program  has  been  written  and  tested  which  allows 
the  computer  to  assume  the  role  of  both  players.  This  means,  for  example, 
that  the  "machine  learner"  can  play  a  number  of  games  versus  a  machine  opponent 
having  a  fixed  level  of  skill.  Or  that  both  machine  players  can  learn  the 
game  but  use  different  learning  techniques.  The  program  is  constructed  to 
facilitate  the  substitution  of  various  learning  techniques.  Through  a  pecul¬ 
iarity  of  the  game  it  is  possible  at  any  point  in  the  learning  process  to 
determine  objectively  the  level  of  skill  of  a  machine  player,  thus  allowing 
various  learning  techniques  to  be  compared  easily. 
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3.4.2  Electronic  Magic 

In  recent  months  computer  ^Deration  of  musical  sounds  has  ceased 
to  be  merely  a  novelty  and  has  become  a  serious  research  effort  (at,  among 
other  places,  Bell  Telephone  Laboratories).  There  are  many  reasons  for  this, 
among  them  the  fact  that  computers  can  generate  sounds  difficult  or  impossible 
for  human  musicians  to  create  and  can  sight  read  perfectly.  This  latter 
feature  allows  a  composer  to  try  complex  musical  ideas  without  allowing  a  long 
rehearsal  time  for  musicians. 

The  multiple  accumulator  design  of  the  CSX-1  computer  lends  Itself  very 
well  to  the  real-time  generation  of  multi-voice  music.  Circuits  have  now  been 
constructed  which  extend  this  facility  so  that  the  composer  has  control  of 
pitch,  duration  and  loudness  of  individual  notes  (programmable),  attack  and 
timbre  characteristics  of  the  various  voices  (switchable) . 

This  work  was  done  jointly  with  Dr.  L.  A.  Hiller  and  Robert  Baker,  both 
of  the  Music  Department. 

J.  Divilbiss 

3.3  Air  Traffic  Control  Studies 

The  second  of  two  reports  entitled  Automatic  Air  Traffic  Control  will  be 
published  in  March.  Part  II,  CSL  Report  R-l46,  is  sub-titled  "An  Experimental 
Control  Logic."  An  abstract  of  the  second  report  follows. 

The  Work  in  automatic  air  traffic  control  conducted  at  CSL  has  consisted 
of  three  phases:  1)  design  of  an  automatic  control  system,  2)  design  of  an 
experimental  control  philosophy,  and  3)  evaluation  of  the  control  philosophy 
in  real-time  simulation.  This  report,  l.e.  Part  II,  describes  the  experimental 
control  philosophy  and  the  results  of  the  simulation.  The  simulation  was  done 
on  the  Cornfield  Sjetem  which  consists  of  a  special  purpose  tracking  computer 
(TASC)  and  a  general  purpose  digital  computer  (ILLIAC)  connected  together. 

Three  main  control  programs  have  been  written  as  follows. 

1)  PILOT  simulates  manual  control  of  a  square  en  route  area  128  miles  on 
each  side  with  sixteen  airports. 

2)  PICON  (for  PILOT-CONTROL)  simulates  automatic  control  of  en  route, 
approach  and  terminal  sequencing  areas.  En  route  area  is  square,  126  miles  on 
each  side  and  contains  one  airport  with  associated  approach  area  (30  mile 
radius ) . 
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5)  PICCN-ER  (for  PILOT  CONTROL  EN  ROUTE)  simulate*  automatic  control  of 
en  route  and  approach  axe  as.  Control  area  is  square,  512  miles  on  each  side 
with  up  to  sixteen  airports  and  approach  areas. 

In  all  cases,  simulated  traffic  is  made  up  of  k  different  types  of  air¬ 
craft  ranging  in  speed  from  151  aph  cruise  to  5^2  mph  cruise.  The  method  of 
steering  is  radar  vectoring  via  digital  data  link.  Cornfield  System  compu¬ 
tation  speeds  have  limited  syBtem  capacity  to  25  controlled  flights  and  short¬ 
age  of  memory  space  has  limited  the  number  and  complexity  of  the  rules  used 
to  control  traffic.  It  vas  also  necessary  to  assume  zero  winds  aloft,  well 
behaved  air  traffic,  and  one  altitude.  Hie  Cornfield  System  limitations  are 
shown  to  be  non-existent  in  the  system  described  in  Part  I. 

In  the  automatic  control  programs,  the  control  philosophy  in  en  route 
phase  of  flight  consists  of  navigation  and  conflict  resolution.  Navigation 
consists  of  steering  traffic  directly  toward  respective  destinations  and 
conflict  resolution  consists  of  steering  to  avoid  conflicts  by  a  defined  margin. 
In  extreme  cases,  conflicts  are  resolved  by  resorting  to  a  standard  1  minute 
holding  pattern.  In  approach  area,  navigation  is  the  same  as  described  above 
but  the  only  method  of  resolving  conflicts  is  by  holding.  In  terminal  sequenc¬ 
ing,  a  passive  form  of  conflict  avoidance  is  employed  and  traffic  is  steered 
only  to  provide  defined  spacings  in  landing. 


L.  Kypta 
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4.  SYSTEM  THEORY 
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4.1  Introduction 

Ibis  group  Is  concerned  with  three  general  areas  of  research.  Che  area 
Is  automatic  control  theory.  For  the  most  part,  attention  Is  focused  on 
various  aspects  of  optimum  and  adaptive  control  systems.  Another  area  of 
research  Is  circuit  theory  Including  linear  graph  theory  and  computer  methods 
In  network  design.  Applications  of  linear  graph  theory  to  communication  nets 
and  switching  theory  are  being  studied.  A  third  area  of  research  Is  sequential 
machine  theory.  Presently,  work  Is  In  progress  on  algorithmic  design  of  sequen¬ 
tial  switching  circuits,  automata  theory,  and  self -repair. 

4.2  A  Method  of  Control  Law  Generation  Using  Sequential  Phase  Space  Constraints 
Progress  during  the  last  quarter  may  be  broken  up  Into  three  parts.  These 

parts  are  listed  below  with  a  brief  description  of  the  work  done. 

4.2.1  Generalization  of  Constraint  Mapping  Technique 

12  3 

As  reported  earlier  ’  '  a  constraint  mapping  technique  has  been 
devised  for  mapping  state  space  constraints  Into  the  control  space  and  taking 
the  intersection  of  this  subset  with  the  control  space  subset  (given  as  part  of 
the  problem  statement)  to  find  an  admissible  set  of  controls.  It  was  also 

3 

shown  that  the  subspace  so  determined  necessarily  contained  all  admissible 
controls.  However,  It  was  not  sufficient  to  use  a  control  from  this  subspace. 
Further  it  was  shown  that  one  could  define  a  region  of  the  control  space  which 
did  contain  exactly  all  admissible  controls  but  that  this  region  was  a  function 
of  the  present  state  of  the  system.  After  further  study  it  turned  out  that 
this  procedure  could  be  applied  to  a  much  wider  class  of  systems  than  originally 
thought.  The  class  to  which  the  technique  applies  is  called  the  class  of 


Quarterly  Progress  Report,  June,  July,  Aug.,  1962. 

Quarterly  Progress  Report,  Sept.,  Oct.,  Nov.,  1962. 

R-155,  "A  Constraint  Mapping  Technique  for  System  Optimization, "  Nov.  5, 
1962. 
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Invertible  systems.  Two  subclasses  are  defined.  A  careful  definition  of 
this  class  of  systems  follows: 

4.2.2  Definition  of  Invertible  Systems 

If  x(t)  -  f[x(t),  u(t),t]  (1’) 

we  say  that  (l')  Is  Invertible  if  there  exists  an  analytical  expression  fcr  g 
such  that 

u(t)  =  x(t),  t].  (2’) 

Since  we  are  dealing  with  a  problem  of  numerical  solution  of  optimization 
problems,  we  immediately  turn  our  attention  to  the  discrete  version  of  (1*) 
and  (2 1 )  namely 

x(K4*)  =  f[x(K),  u(K),  K,  A]  (1) 

and 

u(K)  =  ffx(K-tA),  x(K),  K,  A].  (2) 


Two  general  classes  of  invertible  systems  may  be  defined  as  follows. 

A.  Class  a  invertible  systems 

This  class  is  made  up  of  two  subclasses: 
a)  Class  (multi-linear)  invertible  systems: 

x(K-tA)  =  f[x(K),K,A]  +  B[x(K),K,A]  h[u(K)]  (3) 

where 


B[x(K),K,A]  1b  a  n  X  n  matrix  with  the  typical  element 


b1(J  =  b1(J[x(K),K,A], 


h  is  an  n-vector  with  the  typical  element 

n 


auVK)- 


<*> 

(5) 


b)  Class  Og  (non-linear,  non  Interacting)  Invertible  systems: 

x(K+A)  -  ffx(K),K,A]  +  B[x(K),K,A]  H[u(K)], 

where  B  is  as  In  (4)  above  and  H  is  an  n  x  n  diagonal  matrix  with  the  following 
typical  elements: 


lfl. 


huCu^K)]  ,  i«J, 

0  ,  i/j. 


(7) 


B.  Class  B  invertible  systems:  This  class  contains  all  invertible 
systems  not  contained  in  class  ». 

Introduction  of  this  class  of  systems  leads  to  a  new  problem  statement 
for  this  investigation.  The  nev  problem  statement  follows: 

4.2.?  Statement  of  the  Problem 

Given  the  following  class  a  invertible  system 


*(t)  *  f[x(t),t]  +  B[x(t),t]  hfu(t)] 

determine  u(t)  such  that 

F[x(t),  u(t),t]  is  minimized,  and 
*"(t)  <  x(t)  <  x+(t),  and 
u"(t)  <  u(t)  <  u+(t). 

The  terms  are  identified  as  follows: 

x-(t) 
u(t) 


(n-dimensional)  state  vector, 
(n-dimenslonal)  state  space  constraints, 
(m-dimensional)  control  vector. 


(8) 

(9) 

(10) 
(n) 


u-(t)  (m-dlmenslonal)  control  space  constraints, 

F  performance  index, 

f  (n-dlmenslonal)  system  vector, 

B  (n  x  n  dimensional)  interaction  matrix. 

The  control  vector  u°(t)  which  satisfies  (9)  subject  to  (8),  (10),  (11), 
is  called  the  globally  optimum  control  vector.  If  u'(t)  is  such,  for  all 
6u(t)  a  small  perturbation  in  u(t),  that 

F[x'(t),  u'(t),t]  <  F[x ' (t )  +  6x(t),  u'(t)  +  6u(t),t], 

then  u'(t)  is  a  locally  optimum  control  vector. 
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The  techniques  of  this  study  determine  a  locally  optimum  control  vector. 

To  find  the  globally  optimum  control  vector,  we  simply  select  the  "best" 

12  K 

locally  optimum  one.  By  this  Is  meant  that  If  u(t )  u(t ),...,  ji(t)  are 
locally  optimal,  then  if  there  exists  J  such  that,  for  all  i/j, 

F[Jx(t),  ^u(t),t]  <  F[^x(t),  ^(tjjt]. 

We  call  ^u(t)  the  globally  optimal  control  vector.  Since  the  number  of 
locally  optimal  control  vectors  may  be  very  large,  generally,  as  In  this  work, 
one  is  satisfied  to  find  a  locally  optimal  control  vector. 

4.2.4  Computer  Results 

The  computer  programs  which  mechanize  the  numerical  procedure 
developed  in  this  Investigation  have  been  written  and  several  examples  run. 

The  results  are  encouraging  and  have  pointed  the  way  to  several  minor  changes 
in  the  optimization  procedure.  The  changes  were  made  and  appear  to  verily  the 
theoretical  predictions  of  behavior  of  the  method. 

4.2.5  Preparation  of  Final  Report 

This  project  will  be .  completed  during  the  coming  quarter .  A  final 
report  has  been  initiated  and  about  one  quarter  of  it  has  been  completely 
written.  The  title  of  this  report  will  be  "On  Direct  Fixed- Time  Optimization  of 
invertible  Systems."  The  fixed-time  designation  Indicates  that  the  control 
period  is  to  be  assumed  fixed  at  the  outset  of  the  optimization  procedure. 

S.  J.  Kahne 

4.3  Linear  Multivariable  Systems  -  Stability  Considerations 

A  method  to  design  optimum,  stable  controllers  was  given  in  the  previous 
progress  report.  During  this  quarter,  the  system  stability  with  respect  to 
plant  parameter  variations  was  considered.  The  necessary  and  sufficient 
conditions  for  the  stability  of  linear  multivariable  systems  are  under  investi¬ 
gation  . 

R.  Narayanasamy 

4.4  Control  Optimization  with  Multirange  Constraints 

The  formulation  of  a  new  control  problem  was  carried  out  during  the  past 
quarter.  We  consider  a  sequence  of  vector  differential  equations  defining  a 
dynamic  process  over  various  ranges  in  time.  The  sequence  of  state  vectors 
may  have  some  components  which  are  common.  We  Impose  control  vector  and  state 
vector  Integrals  constraints  over  the  various  ranges.  The  performance  index 
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nay  Involve  a  sequence  of  state  vectors  and  a  sequence  of  control  vectors.  A 
problem  under  investigation  Is  to  choose  a  sequence  of  control  vectors  which 
satisfies  the  constraints  and  the  sequence  of  process  differential  equations, 
and  which  optimizes  the  performance  Index.  We  have  reformulated  the  problem 
Into  one  which  Involves  multiterminal  constraints  rather  than  Integral 
constraints . 

J.  B.  Cruz,  Jr. 

S.  J.  Kahne 

4.5  Sensitivity  Studies  for  Nonlinear  Feedback  Systems 

This  is  a  new  project  Initiated  during  the  last  quarter.  The  main  objective 
is  to  study  the  effect  of  parameter  variations  In  nonlinear  plants.  The  basic 
question  Is  to  determine  under  what  conditions  feedback  compensation  Is  superior 
to  cascade  or  open  loop  compensation,  from  the  point  of  view  of  plant  parameters 
variation.  Block  diagram  models  for  the  effect  of  plant  parameters  variations 
on  the  complete  Input -output  characteristics  have  been  obtained  using  the  theory 
of  nonlinear  operators. 

J.  B.  Cruz,  Jr. 

W.  R.  Perkins 

4.6  Design  of  Ladder  Filters  by  Computer 

In  19^1,  C.  A.  Desoer  and  others  used  a  steepest  descent  method  to  min¬ 
imize  the  sum  squared  difference  between  the  magnitude  functions  of  Ideal  and 
lossy  networks  for  a  set  of  frequencies.  The  networks  considered  are  passive, 

ladder  filters  with  time  Invariant  elements.  We  plan  to  use  the  Taylor  Series 

5 

approach  which  does  not  deal  directly  with  minimizing  the  norm  of  an  error  but 
minimizes  the  sum  of  the  squares  of  the  difference  between  the  ideal  response 
and  Its  linear  approximation  at  the  sampled  points  (or  discrete  frequencies).  Ex¬ 
perience  has  shown  that  this  approach  Is  powerful  when  the  function  depends 
linearly  on  system  parameters  (e.g.,  network  elements).  If  the  function  Is  not 
linearly  dependent  on  system  parameters,  then  a  modification  such  as  a  relaxation 
method  will  be  needed. 


C.  A.  Desoer  and  S.  K.  Mltra,  "Design  of  Lossy  Filters  by  Digital  Computers, " 
IRE  Trans,  on  Circuit  Theory,  vol.  CT-9,  PP-  192-201;  Sept.,  1961. 

M.  R.  Aaron,  "Tie  Use  of  Least  Squares  In  SyBtem  Design,"  IRE  Trans,  on 
Circuit  Theory,  vol.  CT-3,  pp.  224-231;  Dec.,  1956. 

H.  Sinozakl,  T.  Murata  and  others,  "Application  of  Digital  Computer  to 
Network  Design,"  (In  Japanese),  Convention  Record  of  Japanese  Professional 
Research  Committee  on  Network  Theory,  July,  1962. 


21 


The  method  to  he  studied  requires  the  evaluation  of  the  magnitude  at  a 
set  of  discrete  frequencies.  A  preliminary  version  of  this  method  has  been 
programmed  in  FORTRAN  for  the  CDC  l6o4.  Results  vill  he  reported  later. 

Tadao  Murata 

4.7  Contact  Network  Synthesis 

She  computer  program  for  synthesizing  absolutely  minimal  contact  networks 
for  given  Boolean  functions  is  now  practically  complete  and  is  expected  to  he 
in  production  very  soon.  It  is  hoped  that  the  program  will  provide  answers  to 
some  outstanding  academic  questions. 

W.  Mays da 

S.  Seshu 

4.8  Automata  Theory 

The  literature  survey  undertaken  earlier  in  the  year  is  now  complete  and 
a  report  together  with  an  annotated  bibliography  is  being  Issued. 

H.  Y.  Chang 

4.9  Self-Diagnosis 

A  new  project  was  initiated  during  the  last  quarter  on  the  diagnosis  of 
sequential  machines.  The  general  problem  has  been  discussed  by  Seshu  and 
Freeman*  The  primary  tool  of  this  study  is  a  computer  program.  This  program 
is  currently  being  written  for  the  CDC  l6o4  computer.  This  project  has  four 
major  objectives: 

1.  Understanding  methods  of  programming  a  computer  to  execute  decision 
processes  effectively  and  efficiently  (heuristic  programming) . 

2.  Study  the  problem  of  interconnecting  a  set  of  solutions  obtained  for 
the  pieces  of  a  computer  to  produce  a  set  of  system  diagnostics. 

3>  Study  the  problem  of  organizing  a  digital  computer  to  facilitate  self- 
diagnosis  . 

4.  Provide  the  basic  tool  for  several  future  investigations  (Fh.D. 
theses)  in  this  area. 

The  "working  core"  of  the  program,  comprising  4000  machine  language  ins true 
tions,  is  now  complete.  The  final  program  is  expected  to  be  about  10000  instruc 
tlons  for  the  CDC  l6o4  computer.  Plans  are  under  way  for  encoding  the  logic  of 
the  CSX-1  computer  so  that  the  study  can  utilize  a  "real  system"  as  an  example. 

S.  Seshu 
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4.10  Design  of  Sequential  Switching  Circuits 

Early  computer  results  for  algorithmically  designing  sequential  switch¬ 
ing  circuits  from  input -output  specifications  showed  that  an  unreasonably  large 
number  of  input-output  entries  were  necessary  to  completely  specify  same  siiqple 
problems.  A  more  feasible  specification  appears  in  most  cases  to  be  a  flow 
table.  The  procedures  developed  for  designing  sequential  circuits  from  input- 
output  sequences  were  modified  so  as  to  provide  for  a  flow  table  problem  specif¬ 
ication.  A  machine  language  program  for  the  CDC  l6o4  computer  was  written  to 
verify  the  capabilities  of  the  method.  Results  as  of  this  date  indicate  that 
the  procedure  compares  fairly  well  with  other  procedures,  notably  that  developed 
by  Huffman.  The  number  of  secondary  variables  required  for  realization  varies 
from  1.0  to  about  1.7  times  the  number  needed  for  a  Huffman  realization  in  the 
examples  already  worked,  but  time  estimates  indicate  that  this  procedure  is 
several  orders  of  magnitude  faster  than  that  of  Huffman.  One  example  was  worked 
in  which  varying  quantities  of  the  data  were  available,  the  procedure  gave  a 
realization  requiring  only  1.7?  times  as  many  secondary  variables  as  that  ob¬ 
tained  by  Huffman's  procedure.  (Huffman's  procedure  requires  all  of  the  data 
to  be  available.)  Vith  15  o/o  of  the  data  available,  only  1.25  times  as  many 
secondary  variables  were  needed. 

Other  examples  are  being  prepared  and  some  small  changes  are  being  made  in 
the  program  to  provide  a  better  analysis  of  the  capabilities  of  the  procedure. 

J.  Elsey 

4.11  Self -Repair 

A  program  has  been  written  for  the  CDC  l6o4  computer  to  simulate  various 
organizations  of  systems  with  repair  and  replacement.  It  is  proposed  to  expand 
the  program  to  simulate  self -repairing  systems  composed  of  several  kinds  of 
basic  components.  By  injecting  random  failures,  it  will  then  be  possible  to 
evaluate  the  efficiency  of  the  different  organizations  in  increasing  system 
lifetime  at  the  cost  of  greater  complexity  and  an  increased  number  of  repairs. 

In  its  present  state  of  development,  the  program  is  used  to  simulate  sys¬ 
tems  consisting  of  only  one  kind  of  basic  machine.  Three  or  mom  of  these  work 
in  parallel  with  failures  being  detected  by  a  majority  organ.  When  a  failure 
is  detected,  the  failed  machine  is  replaced  by  a  stand-by  (which  may  also  have 
fslled  while  waiting  to  be  put  into  operation) .  A  time  of  at  least  TQ  is 
required  to  repair  a  failed  machine  and/or  prepare  the  next  stand-by  machine 
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for  operation.  Tc  nay  be  specified  by  an  arbitrary  probability  distribution. 

Hie  results  of  simulating  these  systems  are  being  compared  to  analytical 
calculations  of  system  reliability  In  an  attempt  to  develop  reliable  and  work¬ 
able  exact  methods  of  calculating  the  reliability  or  In  checking  the  validity 
of  any  approximations  made.  In  the  case  of  the  simplest  system  (Tq  »  «*)  a 
solution  has  been  found  which  agrees  excellently  with  the  simulation  results. 

Die  program  has  already  been  found  extremely  valuable  In  Indicating  errors  In 
the  calculations.  It  is  expected  that  its  usefulness  will  even  increase  as 
more  complex  systems  are  studied. 

J.  Kruus 

4.12  Linear  Graph  Theory 

A  digital  computer  program  for  testing  the  realizability  of  a  given  matrix 
as  the  fundamental  cut  set  matrix  of  a  non-orlented  graph  has  Just  been  completed. 
Ciis  is  to  be  combined  with  Seshu's  routine,  whose  purpose  is  to  generate  all 
switching  functions  which  are  equivalent  to  a  given  switching  function.  Bius 
it  will  be  just  a  matter  of  time  until  minimum  contact  networks  are  found  for 
several  switching  functions  whose  minimum  realizations  have  been  unknown  up  to 
the  present. 

W.  Mayeda 

4.15  Communication  Nets 

If  there  exists  one  communication  net  which  has  a  given  matrix  T  as  Its 
terminal  capacity  matrix  then  In  general  there  exist  uncountably  many  such  nets. 
This  Immense  set  of  realizations  of  T  can  be  partitioned  Into  a  finite  number  of 
equivalence  classes.  Among  these  equivalence  classes  for  a  particular  T  there 
may  be  one  which  Is  distinguished  from  the  others  by  virtue  at  having  a  certain 
property.  For  some  time  now  I  have  been  trying  to  prove  that.  In  fact,  this 
special  class  is  present  for  any  realizable  T.  I  recently  found  an  algorithm 
for  obtaining  from  any  realization  of  T  a  member  of  the  special  class  If  it  Is 
nonempty  but,  unfortunately,  this  has  not  yet  shed  any  light  upon  the  general 
question  of  the  existence  of  this  class. 


J.  Beah 
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?•!  Introduction 

The  purpose  of  the  PLATO  project  is  to  develop  an  automatic  teaching  system 
for  tutoring  simultaneously  a  large  number  of  students  in  a  variety  of  subjects. 
The  central  control  element  of  the  teaching  system  is  a  general  purpose  digital 
computer.  The  PLATO  system  differs  from  most  teaching  systems  in  that  a  single 
high  speed  digital  computer  is  used  to  control  all  student  stations.  Thus,  It 
can  bring  to  bear  the  power  of  a  large  digital  computer  in  teaching  each  student. 

5 .2  Physiological  Measurements 

By  the  end  of  January,  electrocardiogram  data  had  been  collected  on  nine 
subjects  vho  worked  completely  through  the  lesson  program  known  as  CHAOS.  Further 
work  is  required  on  equipment  for  recording  physiological  data  during  a  lesson, 
but  this  work  has  been  temporarily  discontinued  in  order  to  concentrate  on  anal¬ 
ysis  of  data  already  acquired,  preparation  of  new  instructional  materials,  and 
a  new  project  (DOPA)  to  be  described  shortly. 

Several  different  modes  of  data  analysis  were  partially  carried  through  in 
order  to  discover  which  mode  of  presentation  maximizes  the  clarity  of  relation¬ 
ships  between  stimulus  and  response  data,  on  the  one  hand,  and  heart  rate  on  the 
other.  It  was  finally  decided  to  represent  heart  rate  as  number  of  beats  per 
tenth  of  a  minute.  Heart  rate  was  then  plotted  as  a  graph  on  the  same  sheet 
as  a  time  response  graph  for  keyboard  data.  One  of  these  graphs  is  shown  in 
Figure  5-1.  Graphs  of  keyboard  data  have  made  it  possible  to  analyze  more 
easily  the  type  of  Intellectual  activity  being  carried  on  by  the  student,  at 
least  for  the  lesson  time  covered  by  keyboard  activity.  The  time  scale  was  div¬ 
ided  into  intervals  representing  discovery,  rediscovery,  false  discovery,  routine 
calculations,  and  indeterminate  activity.  Work  is  now  underway  in  averaging 
heart  rate  for  each  of  these  different  types  of  time  Intervals.  Preliminary 
examination  of  mean  heart  rate  seems  to  show  a  reliable  difference  in  the  level 
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of  activity  between  discovery  periods  and  routine  calculations.  More  definite 
conclusions  await  completion  of  computation  of  variances. 

The  original  CHAOS  program  permits  students  to  discover  patterns  In  num¬ 
erical  sequences.  Work  Is  now  underway  preparing  a  program  which  will  Include 
patterns  in  geometrical  arrays,  such  as  the  triangular  numbers,  and  which  will 
permit  study  of  the  extent  to  which  patterns  so  discovered  can  be  transferred 
to  other  problems  and  also  used  In  combinations  in  a  kind  of  derivation. 
Programming  of  this  material  will  be  carried  out  In  such  a  way  as  to  minimize 
the  periods  in  which  the  students'  activity  is  indeterminate. 

A  new  lesson  program  has  been  launched  to  develop  a  Demonstration  of 
Proofs  in  Algebra  (DOPA)  to  be  performed  by  a  student  In  the  PLATO  lab  mode. 

The  specifications  for  the  DOPA  teacher  have  been  worked  out  in  great  detail 
for  this  program,  pages  have  been  typed  for  photographing  on  slides,  and  the 
constant  list  has  been  prepared.  Preparation  of  the  DOPA  parameters  for  use 
with  PLATO  LAB  is  now  under  way.  In  order  to  detect  the  weaknesses  remaining 
in  this  teaching  logic,  it  is  planned  to  test  this  program  with  University 
High  School  students  who  have  completed  a  unit  in  high  school  algebra  which 
treated  such  proofs.  It  is  also  planned  that  a  program  for  teaching  proofs  of 
this  form  will  be  completed  for  use  with  teachers  attending  the  NSF  Simmer 
Institute  this  summer. 

The  general  idea  of  the  teaching  for  DOPA  is  that  the  computer  will  print 
on  the  blackboard  each  line  of  a  proof  once  a  student  has  supplied  the  prin¬ 
ciple,  the  substitutions,  and  essential  copying  instructions  which  justify 
that  line  as  a  conclusion  from  a  preceding  line.  That  is,  the  machine  supplies 
the  logical  rules  and  refuses  to  act  in  violation  of  them;  the  student  must 
supply  the  strategy  or  heuristics  necessary  to  arrive  at  the  desired  conclu¬ 
sion.  The  major  problem  in  this  project  has  been  to  eliminate  the  possibility 
of  drawing  invalid  conclusions .  A  period  of  testing  this  program  will  be 
needed,  however,  in  order  to  insure  that  no  loopholes  have  been  overlooked. 

It  is  also  expected  that  students  working  on  this  program  will  provide  inter¬ 
esting  data  on  the  psychology  of  heuristic  thinking  and  strategy,  and  that 
physiological  measures  such  as  heart  rate  will  show  interesting  changes  as  the 
student  moves  from  analysis  to  synthesis  in  the  development  of  a  proof. 
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It  is  hoped  that  this  demonstration  program  will  be  a  stimulus  to  other 
writers  to  prepare  a  sequence  of  such  lessons  vhich  vlll  develop  the  ability 
to  create  proofs  by  students  vho  have  not  had  previous  experience  of  this  kind. 

J.  A.  Easley,  Jr. 

5.3  PLATO  Laboratory 

The  PLATO  Laboratory  program  vas  completed  In  December  and  a  demonstration 
lesson  using  applications  of  Archimedes  principle  vas  vrltten.  On  January 
first,  the  ILLIAC  computer  vas  no  longer  available  and  so  the  PLATO  II  equip¬ 
ment  vas  modified  for  use  vith  the  l6o4  computer . 

The  PLATO  Laboratory  program  offers  possibilities  for  "teachers"  different 
from  the  PLATO  II  teacher.  To  demonstrate  its  varied  use,  tvo  "teachers"  have 
been  in  the  process  of  preparation  during  the  last  half  of  this  quarter.  One 
of  these,  DOPA,  is  reported  by  members  of  the  UICSM  project  in  another  part  of 
the  PLATO  progress  report.  Almost  completed  is  a  program  for  Dr.  Richard 
Suchman  and  his  associates  for  use  in  the  scientific  inquiry  study  vhich  he  is 
carrying  on  vith  11  and  12  year  old  students.  The  program  allovs  a  student  to 
see  a  brief  film  demonstrating  the  properties  of  a  bimetal  strip  and  then  to 
follov  his  ovn  path,  via  PLATO  Laboratory  through  questions,  ansvers,  and  sim¬ 
ulated  experiments  until  he  arrives  at  conclusions  about  the  basic  principles 
Involved  in  the  bimetal  strip  experiment .  A  record  of  the  student's  pattern  of 
discovery  is  kept  by  the  computer  and  is  provided  at  the  end  of  each  student's 
lesson.  The  PLATO  Laboratory  program  allovs  progress  back  and  forth  through 
the  lesson  material  via  linked  modes.  The  linking  of  the  modes  and  the  calcula¬ 
tions  performed  in  each  mode  are  specified  by  the  parameters  for  each  lesson. 

The  PLATO  Laboratory  program  vith  its  versatile  linked  mode  construction 
has  suggested  the  possibility  of  a  very  general  master  program  for  PLATO  vhich 
vill  allov  non-computer  trained  persons  to  vrite  their  ovn  "teachers"  for  what¬ 
ever  kinds  of  lessons  they  desire.  The  users  vlll  be  able  to  vrite  the  teachers 
in  a  pseudo  language  vhich  the  computer  vlll  then  be  able  to  read,  compile  and 
assemble  into  the  particular  "teacher"  program.  Work  has  now  been  started  on 
the  general  program  and  ve  plan  to  have  it  completed  in  the  next  quarter. 

D*  L.  Bitzer 

E.  R.  Iyman 
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$.4  PLATO  III  Computer  Program 

In  the  Quarterly  Progress  Report  for  June  1962  a  PLATO  III  program  for 
use  with  the  CDC  l6o4  was  outlined.  During  the  past  quarter  a  substantial 
part  of  this  program  was  coded  and  it  is  expected  that  the  complete  program 
will  be  ready  for  checking  with  the  new  PLATO  III  equipment  by  April,  1963. 

While  coding  of  the  program  was  in  progress,  we  also  considered  what 
lesson  material  would  be  apt  for  use  with  the  new  program,  particularly  for 
demonstrating  its  more  unusual  and  powerful  features.  In  preparation  are  a 
lesson  in  astronomical  measurement  for  grade  school  children,  and  an  exposition 
of  the  elements  of  permutation  groups. 

P.  Braunfeld 
L.  Morgan 

3.3  PLATO  III  System  Equipment 

PLATO  III,  the  system  now  being  developed,  differs  from  PLATO  I  and  PLATO 
II  in  that  a  maximum  of  32  students  (rather  than  only  two  students)  can  be 
taught  concurrently.1  In  addition,  the  PLATO  III  system  will  have  greater 
flexibility  due  to  an  increased  slide  capacity  and  incorporation  of  a  selective 
erase  for  the  blackboard.  A  block  diagram  showing  the  organization  of  the 
PLATO  III  equipment  is  shown  in  Figure  3.2. 

3.3.1  Equipment  Stared  by  all  Students 

Input -Output  Interface.  This  segment  of  the  PLATO  III  equipment 
is  designed  to  act  as  the  buffer  between  the  32  student  stations  and  one 
input -output  channel  of  the  CDC  l6o4  or  the  CSL  CSX-1. 

The  modules  required  for  operation  with  the  CDC-l6o4  are  nearing  completion 
and  should  be  ready  for  testing  as  a  unit  early  in  the  next  quarter.  Most  of 
the  modules  will  plug  into  a  specially  adapted  rack  (which  is  also  nearing 
completion  at  this  time),  while  two  driver  modules  are  being  located  near  the 
computer  to  allow  semi -remote  operation  (the  student  stations  will  be  two  floors 
away  from  the  computers).  Operation  will  start  with  two  student  stations,  with 
additional  stations  added  two  at  a  time  as  completed. 

Slide  Scanner  (electronic  book).  The  present  slide  scanner  In  use 
with  the  PLATO  II  system  has  been  equipped  with  new  photo -multiplier  pre-amps 


1 


Bitzer,  Braunfeld,  Lichtenberger :  PLATO  II:  A  Multiple -Student,  Computer- 
Controlled,  Automated  Teaching  PevlceT  CSL  Report  1-109,  Octo'ber  19&1. 
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which  have  completely  eliminated  a  tendency  for  weak  oscillations  to  appear 
on  certain  selected  slides.  The  present  slide  selector  along  with  a  second 
duplicate  slide  selector  will  be  made  ready  for  use  with  the  PLATO  III  system 
during  the  coming  quarter. 

X-Y  Character,  Sweep,  and  Sweep-sync  Generators.  Prototype  char¬ 
acter  D-A,  TV  sweep,  and  sweep  sync  circuitry  providing  for  write,  read,  blank¬ 
ing,  and  sync  signals  common  for  all  student  stations  is  nearing  conpletion. 

It  is  expected  that  performance  evaluation  and  If  necessary,  further  develop¬ 
ment  of  the  circuitry  will  be  undertaken  in  the  coming  quarter.  In  addition, 
construction  of  prototype  X-Y  circuitry  providing  for  selective  erase  will  be 
undertaken. 

5.5*2  Individual  Student  Equipment 

Storage  Tube  Control  Logic.  Cards  for  one  student  station  have 
been  completed  and  installed  in  the  storage  tube  module.  Testing  has  not  yet 
been  completed.  The  cards  for  another  student  station  have  been  started. 

Storage  Tube  Circuitry  (electronic  blackboard).  Final  layout, 
fabrication,  and  checkout  of  a  prototype,  hybrid-vacuum  tube -transistor  storage 
tube  p  re -amplifier  was  begun.  The  performance  of  the  new  amplifier  appears 
satisfactory. 

Performance  evaluation  of  a  hybrid  vacuum  tube -transistor  deflection 
amplifier  required  for  storage  tube  beam  deflection  was  also  undertaken  this 
quarter.  Early  results  indicate  adequate  recovery  and  stability  characteristics. 

Keysets .  Four  keyboards  have  been  obtained  from  Teletype  Corpora¬ 
tion.  The  circuitry  designed  and  constructed  to  allow  their  operation  with  the 
PLATO  III  system  has  been  completed.  Minor  modifications  to  the  keyboards  them¬ 
selves  must  still  be  made. 

D.  Diamond 
B.  Voth 

5.6  Plasma  Storage  Tube 

The  purpose  of  developing  a  plasma  storage  tube  is  to  find  a  less  expensive 
replacement  for  the  present  storage  tube  system. 

During  the  last  quarter  ejqperiments  have  shown  the  feasibility  of  switch¬ 
ing  the  grid  structure  using  an  R.F.  initiated  plasma  to  establish  a  D.C.  path. 
The  intersection  of  the  horizontal  and  vertical  D.C.  paths  select  the  specified 
spot  to  be  initiated  on  the  "screen"  and  raise  the  potential  there  so  that  a 
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smaller  R.F.  signal  may  initiate  and  subsequently  sustain  this  discharge. 

Ibe  experiment  has  been  performed  at  ten  megacycles  using  neon  discharge 
bulbs  to  slmllate  the  grid  structure. 

A  four  by  four  array  (sixteen  spots)  has  been  constructed  but  difficulties 
In  obtaining  ten  megacycle  transformers  and  neon  discharge  bulbs  with  nearly 
uniform  initiating  potentials  prevented  testing  the  array  In  time  for  this 
report. 

R.  H.  Willson 


6.  VACUUM  INSTRUMENTATION 
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E.  Alpert  Wll  Schuemann 

W.  C.  Prothe  Jose  L.  de  Segovia 

Donald  Lee  Hairy  Tomaschke 

6.1  Sticking  Probability  of  on  Molybdenum 

IXiring  this  quarter  the  sticking  probability  of  Ng  on  molybdenum  vas 
again  measured.  In  previous  eaqperiments  it  vas  observed  that  vhen  Ng  vas 
introduced  into  a  system  the  sticking  probability  vould  increase  vith  coverage 
and  then  fall*  This  unusual  result  could  be  explained  by  a  substitutional 
adsorption  involving  CO.  The  vacuum  system  being  used  to  repeat  the  previous 
experiments  differs  from  the  previous  system  only  in  that  an  omegatron  has 
been  added.  This  work  is  still  in  progress. 

6.2  Anomalous  Behavior  of  Ionization  Gauges 

The  investigation  into  anomalously  high  collector  currents  in  ionization 
vacuum  gauges  mentioned  in  the  last  progress  report  has  been  continued.  We 
are  convinced  that  under  certain  circumstances  the  ionizing  electrons,  vhen 
they  are  collected  on  the  grid  surface,  cause  the  formation  of  ions.  A 
fraction  of  these  surface  ions  are  collected  and  hence  can  give  rise  to  ab¬ 
normally  high  collector  currents.  We  have  definitely  shovn  that  the  double 
collector  gauges  of  Redhead  and  Lee  do  not  give  meaningful  readings  vhen  sur¬ 
face  ionization  currents  are  comparable  to  volume  ionization  currents.  The 
probable  maximum  error  introduced  by  surface  ionization  vhen  a  standard 
Bayard -Alpert  gauge  is  running  at  10  ma  is  not  sufficient  to  reduce  the  use¬ 
fulness  of  the  gauge  for  the  vast  majority  of  applications.  Ibis  vork  is  also 
being  continued. 

6.3  Photo  Suppressor  Gauge  Tests 

We  are  continuing  to  use  the  photo  current  suppressor  gauges  mentioned  in 
four  previous  progress  reports.  Our  typical  vacuum  system  base  pressures 
have  been  shovn  to  be  1  X  10 Torr  vithout  any  cryogenic  pumping.  CXir 
confidence  in  the  gauge  is  increasing  and  it  has  already  proved  to  be  a  val¬ 
uable  research  tool  in  the  investigation  of  the  anomalous  ion  current  men¬ 
tioned  above. 

W.  Schuemann 
Jose  de  Segovia 
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7.  PLASMA  PHYSICS 
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7.1  The  Boltzmann  Equation  for  a  Shock  Wave 

The  quarter  has  been  occupied  both  with  tests  of  the  l6o4-  program  that 
make  use  of  Its  high  speed  and  with  the  development  of  additional  test  programs. 

A  series  of  calculations  was  made  for  a  Mach  number  of  1.3^1  for  samples 
of  two  sizes,  with  and  without  the  moment  correction.  Instability  of  the  type 
found  in  the  ILL LAC  calculations  for  this  Mach  number  have  disappeared  but  the 
fluctuations  still  are  very  large.  The  moment  calculations  near  the  cold 
boundary  sometimes  produced  another  type  of  instability.  The  program  was  there¬ 
fore  modified  so  that  moment  corrections  are  automatically  omitted  when  these 
would  produce  sizeable  distortions  of  the  velocity  distribution  function. 

Sven  though  calculations  of  shock  structure  appeared  to  be  successful  in 
the  previous  quarter  for  the  higher  Mach  number  of  3.17,  two  unsatisfactory 
features  of  the  overall  method  remain:  the  large  fluctuations  for  the  low  Mach 
number  case  and  an  apparent  bias  in  the  calculation  of  the  collision  Integral. 
Bach  difficulty  suggests  the  need  for  making  a  close  study  of  the  calculation 
of  the  collision  integral,  (a-bf),  a  need  that  is  obvious  in  the  case  of  the 
apparently  biased  results,  but  which  requires  a  little  explanation  in  the  caae 
of  the  low  Mach  number  calculations.  For  low  enough  Mach  numbers,  the  Monte 
Carlo  method  may  be  expected  to  fall,  because  the  collision  Integral  must  then 
be  calculated  as  a  small  difference  of  two  quantities,  each  of  which  is  subject 
to  statistical  fluctuations.  (Nordsleck's  version  of  the  Monte  Carlo  method 
cancels  part  of  these  fluctuations.)  Finding  the  lower  limit  of  the  Mach 
number  for  which  the  Nordsleck  calculation  of  shock  structure  is  possible  re¬ 
quires  examination  of  the  accuracy  of  the  calculation  of  each  of  the  two  parts 
of  the  collision  integral.  Note  that  a  bias  in  the  calculation  of  the  collision 
integral  near  the  boundaries  might  disturb  the  calculation  more  for  low  Mach 
e  numbers  than  for  high  Mach  numbers,  for  in  the  later  case  it  may,  in  effect, 
distort  only  the  distance  scale  within  the  shock. 
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Ia  order  that  suitable  tools  to  study  the  collision  Integral  calculation 
would  be  available,  the  second  half  of  the  quarter  vas  devoted  to  vriting  and 
testing  nev  and  modified  routines  (ANLTABF,  DFTAB,  PSIAB,  IOSLET,  NUN,  and  NUP) 
and  incorporating  them  into  the  main  program.  In  various  combinations  these 
new  routines  make  possible  the  following  studies: 

a)  Calculation  of  exact  values  of  a(v)  =  an(v)  for  both  the  cold  and  hot 
equilibrium  distributions; 

b)  Isoline  displays  via  ISOLET  of  values  of  a,  aQ,  a-bf,  a-a^,  f (a/aQ) 
as  functions  of  v; 

c)  Calculation  of  Jx(v)dv  and  Jx(v)v^  dv  where  X(v)  is  alternatively 
a(v)  and  (a-bf). 

These  integrals  represent  such  physically  interesting  quantities  as 
collision  rates  and  density  gradients,  and  in  addition  are  useful  in  monitoring 
calculation  of  the  collision  integral.  The  density  gradient  will  also  be 
necessary  in  the  future  in  an  improved  method  of  integration  of  the  Boltzmann 
equation  itself.  Preliminary  calculations  with  these  two  tools  tends  to  con¬ 
firm  the  bias  in  the  Monte  Carlo  calculation  of  the  collision  integral  and 
suggest  that  possibly  this  bias  is  somewhat  angle  dependent.  At  least  part  of 
the  next  quarter  will  be  devoted  to  searching  for  the  origin  of  these  apparent 
biases.  Monte  Carlo  calculations  of  the  parts  of  the  collision  integral  are 
feasible  on  the  CDC  l6o4  for  40,000  or  more  collisions,  making  possible  val¬ 
uable  studies  of  convergence  as  a  function  of  sample  size. 

A  report  entitled  "Numerical  Studies  of  Strong  Shock  Waves,  Part  II: 
Results  of  ILLIAC  Calculations."  was  prepared  in  October  for  internal  use. 

This  paper,  with  small  modifications,  was  submitted  in  January  for  reproduction 
as  CSL  Report  1-117 . 

J.  K.  Aggarwal 
B.  L.  Hicks 
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7 .2  Momentum  Transfer  Cross  Section  of  Electrons  with  Neon  Atoms 

In  the  preceding  quarter,  measurements  on  the  temperature  (200°  -  600°K) 

dependence  of  the  effective  electron  collision  frequency  in  a  decaying  neon 

plasma  have  been  completed,  from  these,  one  is  able'1'  to  determine  the  energy 

dependence  of  the  momentum  transfer  cross  section  of  electrons  with  neon  atoms. 

2 

Ibis  result,  together  with  the  results  obtained  by  other  workers  at  higher 

energies,  is  presented  in  Fig.  7*1.  Fair  agreement  is  noticed  in  the  overlapping 

range  of  energies,  except  that  the  present  experiment  shows  an  apparent  minimum 

in  cross  section  in  the  neighborhood  of  0.11  electron  volts.  Further  theoretical 

study  on  this  particular  feature  is  proceeding  since  slmost  all  semi -empirical 

3^-5 

theories  *  on  the  scattering  cross  sections  of  slow  electrons  with  neon  atoms 

6  T 

predicted  no  such  minimum.  The  apparent  failures0*  in  observing  a  Ramsauer 

6  T 

effect  by  cross  modulation  *  f  techniques  can  also  be  explained  by  the  present 
experiment  and  analysis.  A  detailed  report  on  this  subject -will  be  prepared  in 
the  near  future. 


C.  L.  Chen 
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Electron  energy  (in  eV) 


Qm,(u)  =  1.315  xlO'16  -1.19  x  10"l5u  +  5.37  x  10‘15u2  cm2 
Qmj(u )  =  1854  x  10' 16  -5.05  x  lCT'V*  +  4.01  x  10'15u2  ci 
Qm,(u)  =  1.07  x  10'16- 2.127  x  10‘lsuS/2+  9.92  x  10' 15  u*2 
where  u  is  electron  energy  in  eV. 


Figure  7.1.  Momentum  Transfer  Cross  Section. 
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1  7»3  Incoherent  Scattering  of  Microwaves  from  a  Plasma 

The  Incoherent  scattering  of  electromagnetic  waves  from  plasmas  Is  receiv¬ 
ing  much  Interest  at  the  present  time  In  conjunction  with  the  backscatter  of 

g 

radio  waves  from  the  ionospheric  F-layer  and  the  possibility  of  using  laser 

beams  for  the  diagnostics  of  high  density  plasmas.  The  theory  of  incoherent 

scattering  of  radiation  from  plasmas  has  been  presented  in  a  number  of 

papers.  *  *  The  characteristic  features  of  the  scattering  cross  section 

are  sketched  schematically  In  Fig.  7.2,  where  the  cross  section  is  plotted 

versus  Ad  =  | to  -  a)Q  | ,  ai  being  the  frequency  of  incidental  radiation.  The 

cross  section  has  a  first  maximum  In  the  neighborhood  of  kv^  (v^  is  the  mean 

Ion  velocity),  a  second  much  smaller  maximum  at  kv  (v  is  the  mean  electron 

velocity)  which  is  emitted  in  the  drawing,  and  a  third  sharp  maximum  at  the 

plasma  frequency.  The  absolute  value  of  the  total  cross  section  is  of  the 
2 

order  r  (r  is  the  classical  electron  radius), 
o  x  o  7 

We  have  made  an  attempt  to  detect  the  incoherent  scattering  of  microwaves 
from  a  laboratory  plasma.  A  schematic  diagram  of  the  apparatus  is  sketched 
in  Fig.  7 *3*  A  microwave  signal  of  the  order  of  50  mw  at  8855  me.  is  prop¬ 
agated  through  the  positive  column  of  a  DC -discharge  housed  in  a  waveguide. 

The  discharge  tube  is  filled  with  neon  at  5  mm  of  Hg  pressure  and  operated  at 
a  current  of  15  ma.  The  electron  density  corresponds  to  a  plasma  frequency  of 
cDp/2jt  ~  1500  me .  The  microwaves  that  have  traversed  the  discharge  are  absorbed 
in  a  termination.  Microwaves  that  are  incoherently  scattered  in  the  backward 
direction  pass  the  circulator  and  are  mixed  with  a  local  oscillator  such  that 
the  difference  frequency  falls  into  the  passband  of  the  IF-amplifier  tuned  to 
50  me,  and  of  about  2  me  bandwidth.  The  microwave  signal  is  modulated  at 
440  c/sec,  by  means  of  a  gyroline,  to  permit  detection  with  a  polarity  coinci¬ 
dence  detector.  Apart  from  incoherent  backscatter ing,  a  considerable  amount  of 
coherent  backscattering  comes  from  the  discharge  tube  and  discontinuities  in 
the  waveguide  system.  A  tuner  is  inserted  before  the  termination  which  permits 
nulling  this  coherent  part  at  the  detector. 
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Pig .  7*4  shows  the  preliminary  results  obtained.  The  output  of  the  pol¬ 
arity  coincidence  detector  is  plotted  versus  the  heat  frequency.  A  slight 
signal  Increase  Is  obtained  In  the  neighborhood  of  +  ojq.  Although  Integra¬ 
tion  times  of  one  hour  were  used,  the  error  is  still  considerable.  More 
accurate  measurements  are  needed  before  a  comparison  with  the  theory  will  be 
attempted.  It  Is  clear  that  the  line  width  at  the  plasma  frequency  will  be 
substantially  Increased  over  the  theoretical  value,  due  to  the  Inhomogeneous 
electron  density  and  the  presence  of  electron -neutral  collisions. 

The  greatest  limitation  In  sensitivity  is  presently  the  noise  of  the 
IF -amplifier.  The  experiment  is  presently  being  modified  In  several  ways  to 
make  more  accurate  measurement  possible. 

An  attempt  to  look  for  the  Doppler  broadening  of  the  central  line  was 
unsuccessful,  due  to  the  large  density  fluctuations  In  dc -discharges  around 
zero  frequency  that  overshadows  the  effect. 

M.  Raether 

W.  Carr  * 

7.4  Plasma  Delay  Line 

The  crossed  electric  and  magnetic  fields  In  the  plasma  delay  line  not 

only  control  the  perpendicular  drift  of  charged  particles,  but  they  also 

constrain  the  particles  from  streaming  along  the  electric  field  lines.  Among 

the  magnetic  lines,  however,  there  is  no  constraint,  and  In  a  practical  device 

some  means  must  be  provided  to  contain  the  particles.  In  the  plasma  mode 

electrons  are  attracted  by  the  ions,  which,  because  of  their  greater  lntertla, 
13 

diffuse  more  slowly  .  In  a  circular  electron  device,  focussing  fields  that 

do  not  vitiate  equal  time -of -flight  conditions  can  be  provided  by  hyperbolic 
Ik 

electrodes.  In  rectangular  geometry,  however,  such  focussing  fields  demand 
the  introduction  of  conductors  In  the  drift  space. 

One  possible  electrode  configuration  is  shown  in  Fig.  7-5,  which  Is  a 
section  of  drift  space  perpendicular  to  the  drift  direction.  The  end  plates 
are  conductors  set  at  fixed  potentials,  and  the  round  dots  represent  cross 
sections  of  conductors  that  terminate  regions  of  uniform  electric  field.  The 
component  of  electric  field,  J^,  which,  together  with  the  magnetic  field 
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Figure  7*^*  Incoherent  Scattering  of 
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Figure  7*5*  Electrode  Configuration  for  Rectangular  Delay  Line. 
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controls  the  drift  velocity,  Is  substantially  uniform  throughout  the  drift 

space.  The  focussing  field,  E  .  on  the  other  hand,  has  two  discontinuities. 

In  the  central  region  It  is  zero*  In  the  adjacent  regions  It  has  a  constant 

magnitude,  but  in  each  case  it  is  directed  to  force  the  electrons  to  the  center. 

An  approximation  to  this  structure  has  been  used  in  the  design  of  a  new 

experimental  model  of  the  plasma  delay  line.  A  photograph,  Fig*  7-6,  shows 

that  each  end  plate  has  been  replaced  by  two  sections  which  are  insulated  from 

each  other.  Each  section  is  connected  to  the  corresponding  section  on  the 

other  site  of  the  device  by  forty-six  turns  of  resistance  wire  wound  on  ceramic 

forms  which  are  then  spring  loaded  to  prevent  slackening  of  the  wires.  The 

-3 

wire  diameter  is  1  X  10  inches  and  the  spacing  between  turns  is  approximately 
-3 

11  X  10  inches.  In  the  right  half  of  the  photograph,  a  coated  glass  rod  is 
visible  in  the  drift  space.  Ibis  resistor  collects  electrons,  but  unlike  a 
metal  conductor  it  does  not  distort  the  electric  field.  At  the  left  is  a 
triode  launcher  with  a  slit  anode  to  define  a  sheet  electron  beam.  Below  the 
triode,  but  visible  only  as  it  emerges  from  the  end  plate,  is  a  second  coated 
glass  rod.  This  resistor  which  serves  as  a  source  of  secondary  electrons  is 
removed  when  the  tube  is  operated  in  the  plasma  mode. 

The  model,  which  is  free  of  certain  field  distortions  present  in  the 
original  tube,  is  assembled  mostly  with  small  nuts,  bolts,  and  ceramic  spacers. 
It  is  thus  possible,  by  changing  parts,  to  use  different  kinds  of  collectors, 
and  to  vary  mechanical  parameters  such  as  the  slit  width  in  the  triode.  The 
tube  has  been  mounted  in  the  magnet,  and  the  assembly  of  instrumentation  for 
the  experiments  is  almost  complete. 


H.  G.  Slottow 
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8.  HIGH  MAGNETIC  FIELD  SUPERCONDUCTORS 


C.  B.  Satterthwaite  T.  C.  Leung 

J.  0.  Kopplin  R.  Ries 

R.  Klingbiel 


8.1  Introduction 

This  group  is  primarily  interested  in  studies  of  superconductivity  pertain¬ 
ing  to  the  persistance  of  superconductivity  in  very  high  magnetic  fields.  There 
are  three  phases  of  this  program  in  which  work  is  being  initiated.  They  are: 

(a)  basic  studies  of  superconducting  phenomena,  (b)  preparation  and  study  of 
materials,  and  (c)  devices  and  apparatus  involving  superconductors. 

8 .2  Vanadium  Study 

From  the  single  crystal  of  vanadium,  the  preparation  of  which  was  reviewed 
in  the  last  progress  report,  two  parallel  sections  2.5  inches  in  length  were  cut 
using  the  spark-cutting  technique.  One  section  was  cut  to  a  0.020  inch  diameter; 
the  other  to  a  0.100  inch  diameter.  These  sections  were  cut  by  using  for  the 
opposing  electrode  in  the  spark-cutting  procedure  a  0.25  inch  thick  block  into 
which  a  series  of  holes  of  decreasing  diameter  were  machined.  The  single  vanad¬ 
ium  crystal  was  "passed"  through  the  holes  until  the  desired  diameter  was 
reached. 

A  suitable  procedure  for  working  and  reducing  pure  vanadium  to  a  0.020  inch 
diameter  wire  has  been  developed  by  Fansteel  Metallurgical  Corporation  under  a 
joint  research  agreement. 

Preliminary  cryogenic  tests  of  the  0.020  inch  diameter  single  crystal 
indicate  nearly  ideal  characteristics.  Careful  measurements  of  the  transitions 
temperature,  the  critical  field  and  residual  resistance  are  planned. 

Critical  field  measurements  on  highly  worked  0.020  inch  diameter  specimens 
of  pure  vanadium  indicate  a  highly  irregular  pattern. 

J.  0.  Kopplin 

8.3  Niobium-Tin  Wire 

In  cooperation  with  the  Fansteel  Metallurgical  Corporation  some  studies 
have  been  made  of  means  of  fabricating  Nb^Sn  •wire  for  use  in  superconducting 
solenoids  and  other  apparatus.  The  metallurgical  processing  and  fabrication 
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have  been  carried  out  by  Fansteel,  and  the  evaluation  of  superconducting 
properties  has  been  done  by  the  cryogenics  group. 

Two  batches  of  wire  have  been  fabricated.  A. control  batch,  similar 
to  the  vlre  first  developed,  vas  made  by  packing  a  mixture  of  Nb  and  Sn 
powders  into  a  Nb  tube  and  subsequently  drawing  down  to  0.015  inches  diameter. 
The  other  batch  was  made  by  inserting  into  a  l/4  inch  Nb  tube  a  cable  of  con¬ 
tinuous  wires  of  Nb  and  Sn  and  drawing  this  composite  down  to  0.015  inches. 

It  was  hoped  that  the  cable  technique  would  lead  to  a  large  Nb-Sn  interface, 
continuous  along  the  wire  where  reaction  between  the  two  components  could 
occur  and  that  this  structure  would  carry  higher  superconducting  currents  and 
would  have  more  uniform  and  predictable  superconducting  properties. 

Samples  of  both  batches  of  wire  were  heat  treated  in  vacuum  at  several 
temperatures  between  700°C  and  1,050°C  and  at  times  varying  from  4  to  l6  hours 
These  were  evaluated  by  measuring  the  critical  current  as  a  function  of  mag¬ 
netic  field  in  fields  varying  from  0  to  22  kilogauss. 

Critical  currents  for  the  wire  made  from  the  cable  varied  with  magnetic 
field  approximately  as  expected  for  Nb^Sn  but  were  lower  by  a  factor  of  2  or 
3  than  comparably  treated  control  specimens.  Maximum  currents  in  the  cable 
wire  at  22  KG  were  about  15  amps.  Of  the  time  and  temperatures  selected  for 
heat  treatment  900°C  for  l6  hours  seemed  to  be  optimum. 

Recently  developed  metalographic  techniques  applied  to  these  wires  in¬ 
dicate  that  far  less  intimate  contact  between  Nb  and  Sn  was  achieved  in  the 
cable  than  expected.  Means  of  increasing  the  Nb-Sn  Interface  are  being  sought 

R.  Kllngblel 

R.  Rles 

C.  B.  Satterthwalte 

8.4  Preparation  of  Crystalline  Nb^Sn 

A  number  of  experiments  have  been  conducted  using  the  Induction  furnace 
described  earlier  for  the  preparation  of  crystalline  Nbj  Sn.  In  a  graphite 
crucible  a  melt  of  Sn  was  maintained  at  a  temperature  of  1000°C  to  1250°C  in 
intimate  contacts  with  Nb  granules  for  an  extended  period  of  time  to  obtain  a 
saturated  solution  of  Nb  in  Sn.  A  probe  of  Nb  or  a  quartz  probe  carrying  a 
seed  crystal  of  Nbj  Sn  was  lowered  until  it  touched  the  surface  of  the  melt. 


*7. 

With  the  Nb  probe,  because  of  the  metallic  thermal  conductivity  of  the 
Nb,  a  cool  spot  was  presented  at  the  surface  of  the  liquid  to  provide  a  point 
of  nucleation.  In  all  cases  nucleation  occurred  and  clusters  of  crystal 
formed  -with  some  dendritic  growth  along  the  surface  indicating  that  radiation 
cooled  the  surface.  With  the  seed  crystal  growth  was  still  in  clusters  with¬ 
out  growth  of  any  large  single  crystal,  however,  by  both  methods  crystal  faces 
several  tenths  of  a  millimeter  have  formed. 

X-ray  analysis  of  crystals  grown  from  melts  at  1000°C  to  1250°C  has  in¬ 
dicated  that  the  polycrystalline  clusters  are  pure  Nb^  Sn. 

Work  is  in  progress  to  vary  conditions  to  encourage  growth  of  larger  and 
more  perfect  crystals. 


W.  C.  Leung 


9.  VACUUM  EU2CTRICAL  BREAKDOWN 


E.  M.  Lyman 

F.  Konrad 
T.  Casale 

H.  E.  Tomaschke 


9«1  Refinements  of  Preliminary  Measurements 

The  preliminary  measurements^  on  vacuum  breakdown  between  parallel  tungsten 
electrodes  have  been  repeated  with  great  care  to  Include  the  effect  of  the 
power  supply  ripple  on  peak  voltages  applied  to  the  electrodes  and  the  resulting 
peak  field  emission  currents.  The  Fowler -Nordheim  plots  of  log  i/V2  vs.  1/V 
are,  as  before,  essentially  straight  lines  over  the  entire  range  of  emission 
current  (5  to  6  orders  of  magnitude).  Defining  the  breakdown  voltage  as  the 
applied  potential  difference  at  which  the  first  visible  bright  spark  or  arc 
occurs  between  the  electrodes,  the  critical  field  at  the  cathode  for  each  of  the 
breakdown  points  was  calculated  from  the  breakdown  voltage  and  the  field  enhance¬ 
ment  factor  due  to  surface  irregularities  (as  calculated  from  the  slope  of  the 
corresponding  Fowler -Nordheim  plot).  The  critical  field  calculated  in  this 
manner  was  found  to  be  essentially  independent  of  electrode  spacing  over  the 
entire  range  studied;  5  X  10  J  <  d  <  0.4l  cms.,  with  a  mean  value  of 

6.5  X  10  volts/cm.  The  mean  of  the  previously  reported  values  was 

7.6  x  10^  volts/cm. 

In  some  of  the  cases  there  is  a  deviation  of  the  observed  value  of  l/V2 

from  the  straight  line  predicted  by  the  Fowler -Nordheim  equation  in  the  region 

of  large  currents  Just  prior  to  breakdown.  The  deviation  is  toward  higher 

2 

currents  than  prediced  and  is  in  the  general  direction  observed  by  Goftaan  • 

Cae  greatest  deviation  observed  was  by  a  factor  of  2.  About  half  of  the  F-N 
plots  showed  no  deviation. 

Sharp,  faint  spots  of  blue  light  are  observed  on  the  anode  at  voltages  less 

than  the  breakdown  voltages.  These  are  thought  to  be  due  to  impinging  electron 

beams  which  originate  at  the  surface  irregularities  on  the  cathode.  At  voltages 

2 

Just  below  breakdown  there  may  be  dozens  of  such  spots  per  cm  •  Similar  spots 


Coordinated  Science  Laboratory.  Progress  Report  for  June,  July,  August, 
1962.  p.  72 

I.  I.  Gofman,  0.  D.  Protopopov  and  G.  N.  Shuppe,  Soviet  Physics  Solid 
State  Physics  2,  1203  (i960). 


3 

were  observed  by  Little,  Smith,  and  Arnett  ,  and  point  sources  of  x-rays  coming 

4 

from  extended  anodes  and  cathodes  have  been  observed  by  Singer  and  Doolittle. 

It  seems  reasonable  to  conclude  that  for  extended  electrodes,  the  prebreakdown 
current  originates  from  many  points  on  the  cathode,  and  that  at  breakdown  arc¬ 
ing  occurs  from  one  or  more  of  these  points,  possibly  due  to  the  mechanism 

5  6 

postulated  by  Dyke  and  extended  by  Alpert  and  Lee. 

E.  M.  Iyman 

F.  Konrad 
T.  Casale 

9.2  Electron  Microscope  Studies 

Continued  studies  of  the  cathode  surface  with  an  electron  microscope  have 
shown  the  following: 

1)  There  are  projections  or  "whiskers"  on  the  cathode  surface  (see  Fig.  9*1) • 

2)  One  or  more  of  these  whiskers  disappear  after  electrical  breakdown. 

It  should  be  noted  that  for  these  studies  "breakdown"  is  defined  as  the  point 
at  which  the  field -emission  current  suddenly  changes,  either  increasing  or  de¬ 
creasing,  while  the  voltage  is  constant.  This  change  in  current  is  usually  50 
per  cent  or  greater.  Typically  this  is  not  accompanied  by  a  visible  discharge 
as  in  the  usual  definition  of  breakdown  (as,  for  example,  in  Lyman's  work). 

When  the  cathode  tip  has  been  smoothed  by  heating  to  about  2000°  C.,  the 

field  enhancement  factor  as  calculated  from  a  Fowler -Nordhelm  plot^  is  about 

20.  After  one  or  more  breakdowns  this  factor  Increases  to  50  or  greater  with 

-9 

some  values  as  high  as  500.  The  calculated  emitting  areas  range  from  10  to 
-l4  2 

10  cm  .  Enhancement  factors  and  areas  calculated  from  electron  micrographs 
are  in  rough  agreement  with  those  calculated  from  the  Fowler -Nordhelm  plots. 

Since  an  electron  micrograph  gives  only  one  profile  of  the  cathode  surface  it 
is  possible  that  only  part  of  a  whisker  may  be  seen.  Thus,  accurate  calcula¬ 
tions  cannot  be  made  by  this  method. 


5  Little,  Smith  and  Arnett,  NRL  Report  No.  5671,  Oct.  2,  1961. 

4 

B.  Singer  and  H.  D.  Doolittle.  Private  communication. 

^  W.  Dyke.  See,  for  example,  Fhys.  Rev.  §1 ,  1043  (1953) • 

^  D.  Alpert  and  D.  Lee.  Coordinated  Science  Laboratory  Report  R-129, 

June  7,  1962. 
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There  are  several  factors  or  phenomena  which  are  not  understood  at 
present.  One  of  these  Is  the  large  variation  in  the  current  density  calculated 

at  the  breakdown  points.  The  values  of  current  density  range  from  1  x  10°  to 

0  2  T 

2  x  10  amp/cm  •  According  to  the  mechanism  proposed  by  Dyke  ,  breakdown 

8  2 

should  occur  at  a  critical  current  density  of  approximately  10  amp/ cm  .  The 
large  variation  can  be  partially  explained  by  the  difference  in  the  geometry 
of  Dyke's  model  and  that  of  the  observed  whiskers. 

Another  phenomenon  is  a  break  in  the  Fowler-Nordheim  plot  such  that  the 
plot  consists  of  two  intersecting  straight  lines  as  shown  in  Figure  9.2.  Ibis 
phenomenon  is  found  only  occasionally;  but  when  it  is  found,  it  is  reproducible 
provided  no  breakdown  occurs. 

H.  £.  Tomaschke 
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W.  P.  Byke  and  J.  K.  Troian,  Fhys.  Rev.  799  (1953). 


Figure  9*1  (»)•  Electron  micrograph  of  entire 
cathode  tip. 


Enlarged  view 
of  one  section. 


10.  PHYSICAL  ELECTRONICS 
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R.  N.  Peacock  George  Riddle 

10.1  Evaporated  Thin  Films 

The  principal  effort  during  this  quarter  has  been  directed  toward  study 
of  the  voltage -current  characteristics  of  metal -oxide -metal  film  sandwiches. 

In  contrast  to  the  reports  of  same  investigators,  negative  resistance  has  been 
observed  to  frequencies  of  at  least  several  hundred  cycles  per  second,  and 
with  the  samples  in  air  as  well  as  in  vacuum.  Measurements  have  been  made  with 
SnOg  as  well  as  with  the  Al^O^  used  earlier.  Sandwich  structures  made  with 
insulating  tin  oxide  have  properties  similar  to  those  with  aluminum  oxide.  Tin 
oxide  films  prepared  by  sputtering  tin  in  tank  helium  or  air  have  very  good 
dielectric  properties. 

The  receipt  of  a  Varian  VI-4  vacuum  system  in  March,  1963  is  expected  to 

aid'  in  the  clean  preparation  and  measurement  of  films. 

R.  N.  Peacock 
G.  Riddle 


11.  ARMS  CONTROL  AND  DISARMAMENT  RESEARCH 


B.  L.  Hicks  C.  B.  Satterthwaite 

The  part-time  activities  of  several  members  of  the  staff  during  the  past 
quarter  have  reflected  their  interest  in  encouraging  the  growth  of  ACD  research 
at  the  University  of  Illinois ,  Among  these  activities  were: 

a)  attendance  of  Drs.  Satterthwaite  and  Hicks  at  the  International  Arms 
Control  Symposium  at  the  University  of  Michigan,  December  17-20,  1962; 

b)  formation,  with  the  help  of  faculty  from  several  other  departments, 

of  a  weekly,  interdepartmental  Seminar  on  Arms  Control  and  Disarmament  Research 

c)  discussion  with  scientists  of  the  U.S.  Arms  Control  and  Disarmament 
Agency  about  that  Agency* s  programs  of  research  and  of  research  support. 

One  purpose  of  these  activities  is  to  provide  background  for  finding  ACD 
research  problems  which  conform  to  our  professional  interests  and  capabilities* 

B.  L*  Hicks 

C.  B.  Satterthwaite 
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U.l.  Naval  Elactronlca  Laboratory 
Ian  Diego  J2 .  California 
Attn;  Coda  2(00.  C.R.  Manning 


Commanding  Off  leaf  and  Director 
U.l.  Navy  Electronic!  Laboratory 
Ian  Diego  SI.  California 
Attn;  Library 


I  Raaearch  Anatyala  Corporation 

IDS  Arlington  Road 
Betheede  14.  Maryland 
Awn:  library 

5  Scientific  A  Technical  Information  Facility 

F.  O.  box  S700 
betheeda,  Maryland 
Attn:  NASA  A* pr* .amative  IS -AX/DU 


Office  of  Naval  Reaeerch  branch  OFfloa 

1000  Oaary  Mreet 

Ian  Franclaoo.  California 

btanfotd  Elactronlca  Laboratortaa 
IMnford  Uni  vanity 
Itanfcrd.  California 
Attn:  EEL  Documanta  Librarian 

AFRIT  -  SC 
Heedguertera .  U1AT 
waahlngton  2s,  D,  C. 

Director  of  Science  and  Technology 
Headquarter.,  URAF 
waahlngton  25.  D.  C. 

Attn:  AFRIT-II/DU 


Commander 

Air  Force  Cambridge  Raaearch  Laboratortaa 
Laurence  0.  Hanioom  Field 
Sedford ,  Maaaachuaatta 
Ann:  CRXL 

Raaearch  Laboratory  of  Elactronlca 
Maaaachuaatta  Inatttuta  of  Technology 
Cambridge  2*.  Maaaachuaatta 
Attn:  Document  Room.  21-22? 

Lincoln  Laboratory 

Maaaachuaatta  Inatttuta  of  Technology 
F.O.  box  72 

Lexington  72.  Maaaachuaatta 
Attn;  library,  A-0II 


Headquarter.,  RAT  Divtilon 
•oiling  Mr  Toroe  Sara 
Waahlngton  25.  □.  C. 

Attn;  RTHR 

Heedguertera,  U.l.  Army  Material  Command 
Raaearch  Dtvialon,  RAD  Directorate 
Waahlngton  IS,  D.  C. 

Attn:  Phyalca  A  Elactronlca  Iranoh 
Elactronlca  lection 


Commanding  Officer 

Diamond  Odnenoe  Futa  Laboratortaa 

Waehtngton  Is,  D,  C. 

Attn:  librarian.  Room  111,  lidg.  fl 

Operation!  Evaluation  Ooup 
Office  of  tha  CKO  (OpOlED) 

Navy  Department 
Wbehlngton  21,  D.C. 

Chief  of  Naval  Opera  dona 

Tech.  Anaiyau  A  Advisory  Group  fOb-07T1 

Waahlngton  21.  D,  C. 

Commanding  Offlcor  A  Director 
David  w.  Taylor  Modal  beam 

WaaMngton  7,  D.  C. 

Attn:  Cede  141,  library 

•weauof  ftlpg 
Department  of  the  Navy 
RHtMngton  ll.  D.  C. 

Attn:  Ceda  lb* 

Bweeu  of  Mile. 

Navy  Department 
WtaAtngtan  IS.  D.C. 


1  bylvanla  Electric  Product!  Inc. 

Electronic  by  item  a 
Waltham  Laba.  Ubrary 
100  First  Avanua 
Waltham  St,  Maaaachuaatta 

I  Minima rolli- Honeywell  Regulator  Co. 

Aaronauttoal  Dtvialon 
2SOO  Ridgeway  Aoad 
Minnie  poll.  12,  Minnesota 
Attn:  Mr.  D.  F.  Elwell 
Main  button:  SIS 

1  Importer  of  Navel  Malarial 

bureau  of  Ihlpa  Technical  Ragraaentatf 

1*01  Waat  Minnehaha  Avanua 
•t.  Paul  4,  Mlnnaaou 


20  Activity  Supply  Oflioet 

building  2504 
Charles  Wood  Meg 
Tort  Monatouth.  New  Jaraey 
Marked:  For  Office  of  Engineering  Opera tlona 
Order  He.  m?t-PM-t»-tl-ll 

I  Radio  Oarpentton  of  America 

RCA  Lahoratortaa 
David  lamoff  Raaearch  Center 
brinceton,  New  Mreay 
Attn:  Ubrary 

1  Dr.  j,  H.  Frank 

lil  Wathlngton  Avenue 
Wat  tweed ,  New  Jaraey 


Mr.  A.  A.  bawd i Wow 
■all  Telephone  Uhenlortei 
Recta  IB-127 
WMppeny  Read 
Whlppatsy.  New  Saraey 


I  Mailing  Setae  C-SF 

•all  Aeraayatami  Oempany 
b.  O.  bam  1 
buffalo  S,  New  Vert 
Altai  TaehMoal  library 

I  Cornell  Aeronautical  Laboratory.  Inc. 

44SS  Ooootoo  btraat 
buffalo  11.  Naw  York 
Attn:  I.  b.  Da  amend,  librarian 

I  ■Rdtry  Dytb  sea  pa  Company 

Manna  Dtvtalcii  Ubrary 
1S1  Olen  Cove  Reed 
Carle  Place,  L.l, ,  Naw  York 
Ann;  Mr*,  herhen  Judd 

I  Roma  Air  Development  Center 

Drift  lea  Air  Force  bate,  Naw  York 
Attn:  Documanta  library 
RAALD 

I  Qanarai  Electric  Company 

Advanced  Elactronlca  Cankor 
Cornell  Umveralty 
lthaoa,  Naw  York 
Attn:  F.  Kuahn.  Lttxartan 

1  Ubrary 

light  Military  I  lac  Pen  lea  Department 
Qanarai  Electric  Company 
building  No.  21-101 
•oheneetedy  S.  Now  York 

1  Commanding  Offloet 

U.l.  Army  Raaaaroh  Office  (Durham) 
MU!  CRD-AA-lb,  Mr.  Ulah 
jx  CM,  Dube  Motion 
Durham,  North  Carolina 


I  garmila-DCrtKDCR 

battalia  Mamartel  Inatttuta 
»5  Big  Avanua 
Coltaabua  1.  Ohio 

i  AID  (ASNC) 

Wright -battereon  Alt  rote#  Mae 
Ohio 

1  AID  1AJRH0 

Wtlght- battaraon  Air  Fotoa  boaa 
Ohio 

J  Commanding  Offloar  (AD-SJ 

U.l.  Novel  Mr  Development  Canter 
JolmavlUe.  bmmaylv.nl* 

Attiu  NADC  library 

2  Commending  Oft  tow 
bienbford  Araenal 
Philadelphia  27 .  tarmeylvenle 
Ann:  EM  UFA- 1X0 

1  Donate!  Aoenlca  Corporation 

1071  D*  Haven  Mroet 
trait  Oenahehecken .  Pumnytvuna* 
Attn:  MUa  D.  M.  Raaner 
Librarian 


H.  E.  Cochran 

Oak  Mdga  Nattenel  Laboratory 

b.  O.  ban  x 

Oak  Mdga .  Tannonoo 

Praaldant 

U.l.  Army  AN  Dafanaa  beard 
Fort  Ell...  Tames 

U.l,  AX  Fore*  becurtty  baryte 
Sen  Antonie.  Tenet 
Attn:  OOC-R 

ATriA  Technical  Ubrary  An.  2124 
Arlington  Hall  Mauan 
Arlington  12.  Virginia 
Attn:  TTSU 


Army  Material  Command 
Raaaaroh  DtvleMn 
R  4  D  Dtractarato 
■Mg.  T-T 

Qravelley  taint.  Virginia 


ATMDC  (MDbQb/Capt.  WrlgM) 


